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A "split-field" infrared optical system has been designed and constructed to 
provide simultaneous image pairs in a single frame of an infrared (IR) imager, 
differing only in the direction of linear polarization. The optical train is afocal, 
allowing its use with a variety of infrared imaging devices. The system can operate 
in both long-wave IR ( 8-12 pm) and mid-wave IR (3-5 pm) with interchangeable 
polarizing splitter plates. Previous work at the Naval Academic Center for Infrared 
Technology (NACIT) at the Naval Postgraduate School has demonstrated that 
significant image improvement in infrared image contrast can be obtained by use 
of polarization filtering, especially for targets at sea through suppression of 
polarized sea background emission in the long wave, or of polarized reflection in 
the midwave. That work utilized digital subtraction of sequential image pairs with 
orthogonal polarizations, but suffered from inability to obtain simultaneous 
images and also from problems due to reflections from external polarization filters. 
Both of these problems are eliminated with the new split-field technique. Prelim- 
inary tests of this system with an AGA-780 imager were carried out in a field 
experiment using an Arleigh Burke DDG class ship as target at varying distances 
at sea. For comparison, images were also taken in successive pairs in time with the 
AGA-780 viewing the scene directly ( without the split-field adjunct) with 
interchangeable internal polarizing filters in a rotatable filter wheel. Subjective 
and numerical analysis of the data from the field experiment demonstrated good 
image quality and contrast improvement, and potential for future utility. 
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I. 



INTRODUCTION 



Previous work done at the Naval Postgraduate School has 
shown that infrared image contrast for sea targets can be 
improved by polarization filtering (Refs: 20,21,22). 

Unwanted background radiance which clutters the image can 
consist of emission and reflection components. At grazing 
angles an appreciable part of the water surface emission is 
linearly polarized with the electric vector vertical, while 
direct or indirect solar radiation reflected from the 
surface shows horizontal polarization. Other elements of the 
scene, such as the sky and man-made objects, often emit or 
reflect radiation that is unpolarized or polarized in a 
different direction from that of the sea surface. As a 
result, viewing the scene through a polarizing filter, or 
subtracting two otherwise identical images except for 
orthogonal polarization directions, can improve the target 
contrast against background. 

Polarization filters for the infrared are only recently 
available. These consist of "wire-grid" structures 
deposited by photolithography and etching on transparent 
substrates. Such substrates also need to be antireflection 
coated to reduce unpolarized reflection by the dielectric 
substrate . 
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Special problems arise when dealing with IR as opposed 
to visible radiation. Because all materials emit radiation 
at temperatures above 0 degrees K, sensor background 
radiation often contaminates the final image. Very low 
temperatures must be used at all points from which radiation 
can enter the detector in order to reduce the amount of 
unwanted radiation entering below acceptable levels. Liquid 
Nitrogen and thermo-electrically cooled bodies are commonly 
utilized. Also, the 'narcissus effect' can occur in the 
final image when the cold detector 'sees' a reflection of 
itself from a flat or back-focussing optical surface (Ref. 
23) . Various methods to eliminate these problems are 
implemented in the techniques presented. These are key 
driving factors in the utilization of internal filtering and 
in the design of the Naval Postgraduate School Split-Field 
infrared polarized imaging system. 

Important to the study of image improvement is the 
quantitative characterization of image intensity 
differences. In previous experiments done at the Naval 
Postgraduate School in which the orientation of the 
polarizer had to be changed manually between recording 
orthogonally polarized images, there were changes in the 
image in the time interval during this process 
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(Refs. 1, 21, 22) . Consequently it was impossible to have two 
identical images for comparison. To surmount this problem a 
Split-Image optical system was designed and built at NPS. In 
this system two orthogonally polarized images are recorded 
simultaneously. The images can then be digitized and 
mathematically subtracted to demonstrate the contrast 
improvement. Specific design, implementation, performance 
and image analysis will be explained and presented. 

A description is also given of adaptation of the AGA 
imager with polarizing filters installed in the internal 
filter wheel. This system is useful in recording sequential 
images with orthogonal polarizations. 

This thesis explores the phenomenon of infrared 
polarized radiation, background and target polarization 
characteristics, and experimental techniques for filtering 
and improving image quality. 
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II. BASIC INFRA-RED THEORY 



To fully understand the phenomena associated with 
imaging of polarized infrared radiation it is necessary to 
understand the basic physics used to characterize infrared 
phenomena. This chapter deals with the laws of thermal 
radiation, radiation measurement metrics, and discussion of 
thermal imaging fundamentals. These short presentations 
attempt to give some insight into the selection of the 
infrared band and how it is used. 

A. INFRA-RED TERMINOLOGY AND UNITS 

In the study of radiation, the radiant power exchanged 
between the target, medium, and detector is of interest. A 
set of parameters is defined to describe the radiation at 
each point of interest. Table 2.1 shows some of the 
parameters and their units of measurement (Ref. 19) . The 
last column in Table 2.1 denotes the units in the metric 
system. The entire content of this system of units will not 
be addressed, but some particular units need to be 
emphasized. Temperature is measured in Kelvins (K) . The 
steradian (Sr) unit of solid angle is derived from defining 
the total solid angle around the center of a sphere as 
steradians. Note that these quantities are commonly used in 
inconsistent units. This is done to match imaging system 
dimensions more closely for system performance calculations. 
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Symbol 


Name 


Description 


Units 

(metric) 


M 


Radiant 

Flux 

Density 

(Exitance) 


Radiant Flux leaving 
an element of a 
surface per unit 
area . 


Watts/ cm 2 


E 


Irradiance 


Radiant power per 
unit area incident on 
a surface. 


Watts/cm 2 


I 


Radiant 

Intensity 


Radiant power leaving 
a point source per 
unit solid angle. 


Watts/ Sr 


L 


Radiance 


Radiant power leaving 
or arriving at a 
surface at a point in 
a given direction per 
solid angle per area 
normal to that 
direction. 


Watts/ cm 2 Sr 



Table 2.1 Radiometric Units 

All other physical parameters are measured in conventional 
units. Along the imaging path, M characterizes the power 
leaving the target surface per unit area. Some of this 
radiation is emitted by the target and some is reflected due 
to Irradiance from another source. The distant transport of 
the radiant energy from the target is described by the 
Radiant Intensity , I, which includes the geometrical 
spreading with range. I is the most useful parameter in 
calculating the maximum range of an imaging detector. When 
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the radiation reaches the detector, the Irradiance at the 
receiver of the target radiation is converted into a signal 
for recording. 

B. LAWS OF THERMAL RADIATION 

In order to detect a target, it is necessary to 
characterize the radiation coining from the surface. Most 
infra-red detectors only operate efficiently in a small band 
of wavelengths, so the thermal radiation from a target must 
be known in that range. A set of well known radiation laws 
exists that describe the spectral emittance of thermal 
sources . 

1 . Black Body 

A black body is known simply as an ideal emitter and 
absorber. By Kirchoff's Law, for a black-body the spectral 
absorptivity (a) = emissivity (s) = 1 , so the ref lectance (p) 

and transmissivity (x) must be = 0. For a black body this 

is true at all wavelengths. The spectral radiant emittance 
curve for a black body is shown by Fig. 1, and is described 
by Planck's Law. 
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2 . Planck ' s Law 

Planck's Law is the distribution shown in Fig. 1. It 



is given by 



2nc 2 h/X 5 /rT , 2 

- (W/cm 2 ]iin) 

(e^-l) 



where h is Planck's constant, and k is Boltzman's constant. 
The spectral radiance of the black body can be reproduced by 
the Planck distribution. For non-black body emitters, the 
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radiant exitance is reduced by the emissivity (s) . The 
radiant exitance is then 

^K^hblack-bodv ^ ^ ^ 

with £<1 for all non-black bodies. Examples of emissivities 



Material 


Emissivity 


Polished A1 


.05 (373K) 


Steel Oxidized 


.79 (573K) 


Water ( frost/ ice) 


.98/.96 (2 63K) 



Table 2.2 Emissivities of common materials 

of practical materials are shown in Table 2.2 (Ref. 19) 

If s x (T) = constant independent of wavelength, then the body 

is called a grey body. For most surfaces this is not the 
case, and the spectral emissivity varies greatly with 
wavelength. Many materials behave as grey bodies over a 
limited range of wavelength. 

3. Wien Displacement Law 

In Fig. 1, a maximum radiant exitance occurs at a 
wavelength A^ ax . can be found by differentiating the 

Planck distribution function and equating to zero. The 



9 



result is the Wien Displacement Law. 



X T=2898.18]imK (3) 

Using Wien's Displacement Law, the wavelength with the 
highest spectral radiance is known with the surface 
temperature . 

4. Stefan-Boltzman Law 

If the radiant exitance is integrated over all 
wavelengths for a given temperature, the total power emitted 
per unit area can be found. The equation is found to be, 
for a black body 



M(T)=aT 4 ' (4) 

a is the Stefan-Boltzman constant, 0=5.67 x 1CT 12 W/cm 2 K 4 . 

For other bodies, application of the spectral emissivity 
must be done over all wavelengths. 

C. ELECTROMAGNETIC BANDS AND ATMOSPHERIC WINDOWS 

To understand why imaging is performed in the infra-red 
band, knowledge of the bands and the atmospheric attenuation 
'windows' is necessary. All bodies radiate thermally at all 
wavelengths, so it would appear logical to detect the 
radiation at the peak wavelength defined by Wien's 
Displacement Law. This is not done in every case, due 
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partially to atmospheric effects on the emitted radiation. 

A brief discussion of how the target and atmosphere interact 
follows . 

1. Electromagnetic Radiation Bands 

The infrared region of the electromagnetic spectrum lies 
between the radio (longer wavelength) and the visible 
(shorter wavelength) regions . Fig. 2 shows this region of 
the spectrum (Ref. 19) . Its bandwidth is approximately 
. lpm to 1000|jm. The longest wavelength of interest for 

conventional IR imaging is approximately 15pm since most 
targets and modern quantum detectors operate below this 
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Figure 2. Partial Electromagnetic Spectrum 
wavelength. termed the far infrared limit. The region 
nearest the visible is the near infrared. The "millimeter 
wave" band occurs around 1000pm. 
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2. Atmospheric 'Windows' 

The processes that effect the radiation on its path 
through the atmosphere are extremely complex and will only 
be considered briefly. Attenuation affects the amount of 
radiant power incident on the optical system, which in turn 
affects the overall performance of the system. The 
attenuation mechanisms of primary concern are absorption and 
scattering of the radiation. Limiting or avoiding the 
amount of absorption and scattering in the target path is 
essential for any thermal imaging system because of the 
severity of the effects. Suspended aerosol particles are 
the primary contributors to scattering of the radiation in 
the IR. Aerosols include dust, fog, salt, ash, and many 
other airborne particles. Molecules such as H 2 0 and C0 2 are 
the predominant absorbing species in the atmosphere (Ref. 

19) and define the windows in which IR radiation can pass 
through the atmosphere with any effectiveness. This is due 
to the structure of the molecules and the modes in which 
they can be excited and thus absorb energy. Figure 3 shows 
a plot of atmospheric transmittance versus wavelength over 
the near to far infra-red. This plot was generated by an 
atmospheric propagation modelling code, LOWTRAN 7. LOWTRAN 
7 has resident atmospheric models for aerosol and humidity 
content in addition many other atmospheric characteristics . 
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The model used for the plot is for a Mid-Latitude Summer 
atmosphere, 5km path at sea-level. The 'Rural' designation 
is for the aerosol model describing an air mass over land 
with primary components being dust, sand, and hygroscopic 
particles [19] . 




Figure 3. Atmospheric Windows 



Transmission with no attenuation from the figure has a 
transmittance value of 1. The windows are depicted by the 
regions of maximum transmittance. It is clear that any 
imaging would need to take place for radiation found in 
these regions. Commonly used wavelengths for imaging 
military targets are around 4 . 0|im (short wave) for hot 
sources such as exhaust plumes and power plant components, 
and 8-12|am (long wave) for surfaces nearer to ambient 
temperatures. The relation of the temperature to desired 
wavelength band follows the Planck distribution for maximum 
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radiated power. This correlation between radiated power 
wavelengths and propagation windows must be considered when 
using IR radiation for imaging is desired. 
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III. POLARIZATION PHYSICS 



Polarization occurs in the reflected and emitted IR 
radiation of targets and backgrounds. For purposes of 
enhancing target contrast (as exemplified by ship targets at 
sea) , it is necessary to determine the degree of 
polarization present and its orientation. Fundamental 
physical laws and polarization equations describe the 
important parameters of polarization of the target, and sea 
and sky background. Quantification of these parameters 
leads to expected increases in contrast and is necessary in 
designing and implementing polarizing imaging systems. 

A. STOKES PARAMETERS 

Sir George Gabriel Stokes (1852) discovered that the 
polarization state of any light source could be completely 
described by four observable quantities which have become 
known as the Stokes Parameters. This is done by taking the 
time average of the E field vector amplitude described by 
the polarization ellipse, thus yielding an observable 
intensity (Ref. 17). The resultant four parameters are 

I=E 2 +E 2 

s p 

M=E 2 -E z (5) 

* p 

C=2E E cos5 

S P 

S=2E E sin5 

s p 

E s and E c are the normal and parallel components of the wave 
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electric field, with respect to the plane of incidence. 

For a horizontasl surface the vertical component (p) lies in 
the plane of incidence, also referred to as parallel 
orientation. X denotes the total intensity of the light, M, 
the intensity of linear ’s’ or 'p' polarized light, C, the 
linear +45 degree polarized light, and 5, the circularly 
polarized intensity. 8 is the phase difference between the 

s and p components of the E field. The degree of 
polarization is the ratio of the intensities of the 
polarized to the total intensity [3] . 

1/2 

p= (« 2 +C 2+ 5 2 ) (6) 

X 



The polarization of the radiation to be filtered comes from 
effects due to emission and reflection. Specific laws have 
been developed to describe these phenomena and their 
contribution to the Stokes parameters will be shown. 

B. EMISSION POLARIZATION 

The target and background radiance produced by thermal 
emission from surfaces is inherently polarized [9]. Of 
specific interest for imaging targets at sea is the emission 
polarization of the sea and sky background since the primary 
contrast improvement will be obtained by filtering of the 
background. 
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1 . Sea Emission Polarization 

To find the emission polarization of the sea, the 
surface can be treated as an interface with coordinates 
shown in Figure 4 (Ref. 21) . 



Sun 




Figure 4. Coordinate System of Sea Surface 



Y is the azimuthal direction of the sun. S and D are the 
directions of the sun and observer, respectively. N is the 

wave facet normal, 0 and <|) the zenith and azimuth angles for 
each vector. The light coming from the upper side of the 
interface (+z) consists of a reflected and an emitted 
component, such that the total radiance is given by 
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n(Q)=n p (0)+n s ( 8) 



(7) 



The light emitted from the surface can be considered as 
radiation coming from the sea water side of the interface 
undergoing refraction at the surface. Behavior at the 
interface can be derived from the Fresnel equations and is 
discussed in this fashion by Sandus (Ref. 9) and Sidran 
(Ref. 11) following Reference 16. The emissivity of the 
surface is given by Reference 16. 



2 

e- 

( l + z r ) 



( 8 ) 



z r is dependent on the orientation of the polarization 
(parallel or perpendicular) where: 
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z±— 



s+cos 2 © 



s/2cosQ [ ( s+a-sin 2 ©) 1/2 ] 
s+ (a 2 +i> 2 ) cos 2 © 

z — 

72 COS© [s (a 2 +±> 2 ) + (a-sin 2 ©) (a 2 -b 2 ) + [2ab 2 ) ] 1/2 

a = Re{n 2 ) 

b=Im{n 2 ) 

s=]j (a-sin 2 ©) 2 +b 2 

n is the complex index of refraction of the medium, 0, the 
angle of emission. These equations hold true for any 
specular surface. Dependence on wavelength comes from the 
index of refraction (Ref. 11) . This analysis indicates that 
both the emitted and reflected components of surface 
radiance will show polarization, predominately in orthogonal 
directions. Note that consideration is given only to the 
linearly polarized light. For most materials, there is only 
a small range of angles at which elliptical polarization 
occurs. Even in these cases the phase difference between 
orthogonal components is small enough (due to small 
dielectric constants) that an approximation of linearly 
polarized light at all angles is valid (Ref. 3) . This 
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theory has been applied to the sea surface and proven by 
Sidran (Ref. 11), Baesner and McCoyd (Ref. 8), and Gregoris, 
Yu, Cooper, and Milne (Ref. 20) . . 

In dealing with the sea surface, wave structure alters 
the coordinates of interest in Fig. 4. The wave structure 
is considered to compose a distribution of emitting facets 
with unit direction toward the observer defined by the 
vector s (see Fig. 4), where 



s— s 1 1+ s 2 f+ s 2 k 



( 10 ) 



The orientation of the normals to the facets with respect to 
the normal to the plane defined by s and the z axis is 
defined by the angle a. The effective emissivities of each 
facet referred to the vertical and horizontal directions, 
are given by 



s' =e cos 2 a+£, sin 2 a 

J- J- I 

(ll) 

£',=£ sin 2 a+£,cos 2 a 

i j- i 



The prime symbol denotes the emissivity in the observer 
direction relative to the vertical to the plane containing 
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the z axis and the observer direction s (Ref. 8) . The 
distribution of the slopes of the facets has been 
experimentally determined by Cox and Munk (Ref. 15) . The 
Gram-Charlier equation with empirically derived constants 
adequately describes the distribution of wave slopes (Ref. 
15) . The resulting wave slope distributions are considered 
valid for wind speed <14m/sec and wave slope angle (P)<23 
degrees. Plots of the Cox and Munk distribution are shown 
in Fig. 5 from (Ref. 11) . Because shadowing occurs on the 
back surfaces of the waves in the direction of the observer, 
a correction should be applied to the equations as used by 
Basener and McCoyd (Ref. 8) for percent polarization 
calculation. Saunders (Ref. 10) gives a correction to the 
Basener and McCoyd results. Figure 6 shows the polarization 
of emission for 1 . 0|im wavelength and 14m/s wind speed after 

Sidran (Ref. 11) . The percent of polarization, Q is 
defined by 



0 = 100 * 



£ -£ 
c l x 



£ +£ 
c l x 



( 12 ) 



The results in Fig. 6 hold true for all wavelengths < lO^m 
excluding the absorption bands (Ref. 11). Note that the 
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largest polarization effects are seen at grazing incidence. 

2. Sky Emission Polarization 

The IR radiation from the sky is unpolarized (Ref. 2) . 
Use of polarized filtering components will not contribute to 
the contrast improvement measurements when the sky is viewed 
directly. The sky emission in the image will be subject to 
attenuation by the polarizing filter to the degree of the 
unpolarized transmissivity of the filter. Fig. 7 (Ref. 10) 
shows the sky radiance as a function of zenith angle. 
Specific radiance characteristics can be seen in Fig. 7. 

When viewed toward the zenith, the sky emits less due to the 
fact that the atmosphere consists of water vapor and 
aerosols (absorbers) which also act as good emitters, but is 
less dense. 
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Figure 6. Percent Polarization 
for A,= 1.0 |rm 



Figure 7. Sky Radiance as a 
function of zenith angle from 
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Sky radiance increases with zenith angle toward the horizon 
to that of a black-body with atmospheric temperature. This 
is due to the maximum thickness of the atmosphere 
encountered at that angle and even the small emissivities 
of the absorbers combine (Ref. 10) . 

2. Target Emission Polarization 

Although primary contrast improvement in the at-sea 
scenario will be gained by filtering of sea emission 
polarization, there is polarization due to emission from the 
target itself that must be considered. Painted surfaces 
such as are found on all types of vehicles and ships display 
emission polarization (Ref. 3,4). Figure 8 shows 21 
measurements taken on various paint samples ranging from 
very smooth torough sand-paint mixtures (Ref. 4) . 



— Ideol Surface 




Mtu !• tc « ie ** i 2 * x x x a * jd x w * *a < «c < 



Sample to and Uetre<xemoni Group 

Measured polarized signature components tn the 7.5-12 micron band 
for selected paint samples viewed at 45° from normal. 



Figure 8. Ref [4] 
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The line for an ideal surface in Fig. 8 denotes the amount of 
polarization that would be seen from a perfectly specular 
reflector. The rougher samples show less polarization than 
the smoother. The amount of polarization depends on the 
index of refraction of the paint and any degradation due to 
surface roughness. 

C. REFLECTION POLARIZATION 

For thermal radiation at a surface 

p (0) +a (©) + t ( 6) =1 e ( 0 ) =a ( 0 ) (13) 

As stated previously, the total amount of energy at a 
surface consists completely of the reflected, absorbed, and 
transmitted fractions. For a perfectly opaque body, 
therefore 



e(0)+p(0)=l (14) 

The emissivity (s) was treated in the previous section; the 

reflectance (p) of the sky, sea surface, and target must now 
be considered. Absorption will have an effect on the 
reflectance for each surface considered and will act to 
degrade the reflectance. 
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1 . Sea Reflection Polarization 

The application of the Fresnel Equations to a Cox-Munk 
sea surface slope distribution can be shown to predict 
polarization in reflectance as well as emittance with the 
following relations [ 11] . 



e x =i-P ± e,-i.p, 
£= 1/2 (E^ £, ) 
p= 1/2 (p x+ p, ) 

£ = l-p 



( 15 ) . 



The sea surface acts as an opaque body. Values of s may be 
obtained as discussed in the previous section. The 
percentage polarization seen due to reflection is defined as 
[ 11 ] 



P -P. 
P= 100 * — - — ! 



P A + P, 



( 16 ) . 



2 . Sky Reflectance 

As stated in the emission section, the sky acts as a 
poor emitter near the zenith angle due to the low optical 
depth. Because of the high content of atmospheric water 
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vapor and aerosols, it acts as an excellent absorber. By 
equation (13) it is therefore a poor reflector. 

3 . Target Reflectance 

The characteristics of target reflectance follow the 
same behavior as the wave facet reflectance. The difference 
in the magnitude of the polarized components is a function 
of the index of refraction of the surface material of the 
target. As seen in previous cases, maximum polarization 
effects due to emission, and thus polarization will occur 
nearer to grazing incidence angles. The net polarization 
will be dependent on the balance between thermal emission 
and reflection of incident radiance, related through Eq. 

(15). For the purpose of the experiments performed to 
demonstrate the different polarizing imaging techniques, 
target distance and lack of detailed geometric orientation 
of surfaces precludes exact measurement of target surface 
polarization; however, qualitative effects will be shown. 
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IV. DATA ACQUISITION TECHNIQUES 



To analyze the target radiance as it is affected by 
temperaturer , polarization, and other factors, it is 
necessary to employ thermal imaging equipment for recording. 
The split-field polarimeter and internal filters use the 
AGA-780 thermal imaging system for data acquisition. In 
the case of the internal filters which were designed to fit 
in only the AGA-780 Thermovision, the choice is necessary, 
but the split-field is made to adapt to a number of 
different sensors, particularly the AMBER 4128 focal plane 
array system. All the experimental data taken to be 
analyzed was taken with the AGA-780 Thermovision. 
Accompanying software necessary to digitize and display 
radiance data consisted of the AGA CATS software and 
Interactive Data Language (IDL) programs. The following is 
a brief discussion of the data collection devices. 

A. AGA-780 THERMOVISION 

The AGA-708 Thermovision consists of a two channel 
infrared scanning detector system. The AGA can operate in 
both the 2-5.6(im and 8-12|im bands using liquid nitrogen 
cooled InSb and HgCdTe detectors, respectively. Scanning is 
accomplished by two rotating prisms which scan the image 
across the detector, yielding 1.1 mrad geometric resolution. 
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The input lenses for the Long-wave infrared (LWIR) used for 
the experiments consist of anti-reflection coated Germanium. 
The LWIR band was used for all imaging. The optical system 
is f/1.87, with a 7x7 degree field of view. This f number 
was a key design input for the split system to match the 
aperture of the AGA. The AGA-780 is mounted on a common 
plate with the split-field system with a sliding mount to 
facilitate easy changing between the LWIR and MWIR band 
detectors. Figure 11 in chap. 7 shows the arrangement used 
with the split system. 

B. AGA CATS SOFTWARE PACKAGE 

The analog data provided by the scanning detector is 
digitized for display and measurement by the AGA CATS 
software. The CATS software uses an 8 bit A/D converter to 
convert detector intensity into a 16 color display for 
real-time viewing on a computer screen. Many different 
functions are available to the user such as storing frozen 
image data, image manipulation, temperature readout and 
emmissivity compensation. When using the CATS program to 
provide temperature readings of a target, it is necessary to 
calibrate the system to match target radiance with 
temperature. This is accomplished by recording the thermal 
level seen by the A/D conversion and fitting a curve to the 
data through a least-squares method. The curve is 
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reproduced by the equation [Ref 18] 



ThermalLevel— - OS (17) 

C e B/T -l 

A, B, and C are determined by the curve fitting, OS is the 
offset correction up or down, and T is the temperature 
(ordinate) of the plot. Thermal Level value range 
determines the abscissa. Plots for system calibration for 
the internal filters and with the split-image system added 
are shown in Appendix D. 

The most important feature of CATS is the ability to 
store frozen image data for later analysis. To make 
detailed measurements and perform data manipulation, further 
software is required with greater flexibility. 

C. INTERACTIVE DATA LANGUAGE (IDL) SOFTWARE 

In order to process the image data to a level 
comensurate with the accuracy in measurement desired, a 
software package is required that provides flexibility, 
digital resolution, and data manipulation capability. IDL 
allows the user all of these functions and is required to 
give the data processing power needed. It is produced by 
Research Systems Incorporated, Boulder, CO. Data that is 
stored by the CATS program is read into IDL in array format 
where numerical analysis is performed with the data 
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displayed in an image format. Routines are written to 
perform a number of functions and will be discussed in the 
analysis chapter. 
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V. POLARIZATION FILTERING TECHNIQUES 



Use of polarization signature and background features 
for target discrimination implies the employment of 
polarization filters. The technology for large aperture 
polarizers in the IR is relatively new. Most methods 
commonly used in the visible range are not available in the 
IR. 

A. POLARIZING FILTERS 

Polarizing filters have traditionally been constructed 
using anisotropic crystals or wire grids. Of these methods, 
the wire grid technique has been applied in recent years 
using lithographic techniques. This consists of a 
conducting grid deposited on a transparent substrate. 

Placed in the optical path of the incoming radiation, the 
filter will pass the component of the radiation with the 
electric field orthogonal to the orientation of the grid. 

The grid behaves like a dielectric sheet in this case. The 
component of the radiation with the electric vector parallel 
to the wire grid induces currents in the wires and is 
reflected. Fig. 9 shows a schematic of this process. A 
restrictions on the grid is that the line spacing must be of 
the order of the wavelength of the incident radiation 
(Ref. 6) . Such filters for the infrared are available with 
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clear apertures up to 45mm. Typical transparent substrates 
in the infra-red are materials such as Germanium, Silicon, 
and KRS-5. 




Figure 9. Wire Grid Polarizer 



1. External Filters 

The location of the filter in the optical path has a 
large impact on image attributes. If a flat reflecting 
surface is placed in front of the collimated lens, it is 
possible for the cold detector to 'see' itself in the image. 
This is referred to as the 'Narcissus Effect' (Ref. 23), 
which perturbs the image greatly. Tilting the polarizer may 
alleviate this problem but will reflect radiation into the 
image from outside the field of view. Another drawback to 
using an external filter is the size required to cover the 
aperture defined by the input lenses. It is necessary to 
cover the entire aperture which is a function of the system 
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f number and focal length. Benefits of using an external 
filter are ease of insertion/removal and grid orientation 
changing . 

The large external filter was used in the MAPTIP exper- 
iment in October 1993. The large external polarizer pro- 
vided excellent data for the determination of the degree of 
polarization of the targets and backgrounds and of contrast 
improvement (Ref.l). Conclusions showed percent polarization 
of vertical sea-water emission up to 83% in the long-wave IR 
band. Contrast improvement from emission polarization 
filtering was as high as 36.3%. Large improvements such as 
this justify the utility of the external polarizers and show 
their effectiveness. 

2. Internal Filters 

Filters can also be placed inside the AGA-780 imaging 
system between the main input lenses and the detector. This 
allows hands-off operation and smaller filter size. In 
systems which are configured to use filters in the optical 
path for other types of filtering, small polarizing filters 
can be substituted. These filters will typically be de- 
signed into non-imaging locations. 'Narcissus' spots, 
filter mounting and alignment problems, and other complica- 
tions can thus be avoided with this method. Implementation 
can be accomplished with minor system modification and cost. 
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Installation in an existing filter wheel requires a polar- 
izer for each polarizer orientation to be used, adding to 
system expense. Internal filtering is a highly desirable 
option and experimental use will be discussed in more detail 
later . 

An alternative method of internal polarizing applicable 
to focal plane array systems is placing an independent 
filter closely against each element of the focal plane 
array [6]. Each element of the array can then have any 
preselected orientation of linear polarization built into 
the focal plane mask. An added benefit over rotating 
polarizers is the simultaneous recording of all the selected 
polarizations without manipulating a filter selector. 
Disadvantages are the pixel offset of the image for each 
polarization, and the inherent loss of resolution. The 
practicality of this technique for the 3-5|im band has been 
demonstrated in a NASA- supported program (Ref. 6) . 

B. SPLIT FIELD OPTICS 

A method of polarization in which the polarizer is 
internal to an optical system but external to the imaging 
system itself is the technique of Split-Field polarization 
filtering. The concept of the split field system is that 
orthogonally polarized images can be extracted from the 
f ield-of-view simultaneously with one polarizing filter and 
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presented in adjacent sections of the final image field of 
view. This is accomplished by splitting the beam with a 
polarizer. The Environmental Research Institute of Michigan 
(ERIM) has developed such a system which uses Brewster angle 
reflection to separate one plane polarized component from 
the beam, extracting the orthogonal component from the 
transmitted component through the Brewster plate (Ref. 5) . 
Since the Brewster plate is tilted, a reference is needed to 
keep uncalibrated background from contaminating the image, 
as in the NPS system. Commonality exists between the ERIM 
system and the NPS system. However, the NPS system defines 
the apertures of the two channels and brings them to 
adjacent halves of the final image plane. The ERIM system 
focusses the two channels onto separated areas on a single 
integrated focal plane array. 

The Naval Postgraduate School, Naval Academic Center 
for Infrared Technology (NACIT) has constructed and begun 
testing a split-field system. Polarization is attained 
through a wire grid polarizer tilted to the incoming 
radiation with polarization orientation obstained from grid 
transmission and grid reflection. When considering options 
for selecting an imaging system for the split-field optics, 
the scanning AGA was chosen. The only focal plane array 
available for use at this time operates in the 3-5 (im band. 
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Because of the nature of the targets and environments where 
the split system will be used most, it was decided to first 
build the system to operate in the 8-12 jam region. The 
existing hardware and familiarity using the AGA coupled with 
IDL software and established data analysis techniques made 
it a good choice for examining the new system. The trade to 
go to the AGA is in reduced frame sampling rate and 
resolution. A detailed explanation of the split-field 
technique employed, a system description, and polarizing 
characteristics of the system will follow later. 
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VI. INTERNAL POLARIZING FILTERING WITH THE AGA-780 



Internal filtering of the radiation after it passes 
through the imaging system input lenses can be accomplished 
with the AGA-780 Thermovision Scanning Imaging system by 
using an existing filter wheel. The original purpose of the 
filter wheel was to hold various accessory spectral filters 
in the optical path before the radiation reached the 
detector. Fig. 10 shows the internal layout of the AGA-780. 
One of seven positions of the wheel is selected by rotating 
the knob on the front face of the unit. Care must be taken 




Figure 10. Internal Layout of the AGA 
Thermovision (model with only one IR Band) . 
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to ensure that no internal circuitry compensates for an 
assumed filter present at the position of the polarizer. 
Polarization measurements with this system require 
sequential images with alternate filter wheel positions 
inserting vertically and horizontally polarized filters. 

This necessarily involves time delays. 

A. ALUMINUM/KRS-5 FILTERS 

The filters used consisted of an aluminum grid 
superimposed on a 9.5mm diameter x 3mm thick round KRS-5 
substrate manufactured by Graseby-Specac; Suffolk, England. 
The performance specifications provided with the filters are 
nearly identical for both and will be treated as such for 
explanation. Appendix A shows the performance curves for 
the two filters. The various traces numbered 1-3 represent 
power transmitted with the grid vertical, horizontal, and 
with the grids from two filters crossed, relative to the 
electric vector of the incoming radiation. The determina- 
tion of the perpendicular and parallel transmittances for 
the filter involves use of a spectrometer which itself 
includes internal polarizing reflections. Thus it must be 
accomplished by measuring the transmittance twice, with the 
grid in two orthogonal orientations. These two measured 
transmittances are designated E v and E h . The corresponding 
instrument transmittance' is measured without the grid in 
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place as E 0 , and the transmittance for two crossed grids as 
E x . From these quantities the grid transmittances are 
derived as 



E+E 
K=-± ^ 



* 2 = 



_ 1 



E+E 

h v 



(18) 



These quantities are specified by the manufacturer for each 
filter. The expressions are approximate values derived by 
Young, et al. (Ref. 13). 

B. FILTER INSTALLATION IN THE AGA-780 

Figure 10 shows the internal layout of the AGA-780. To 
install the filters, the appropriate holes in the wheel were 
enlarged approximately 0.5mm in diameter. Two filters were 
then installed in filter positions 4 and 5, carefully 
oriented with passage axis vertical and horizontal, 
respectively. The filter wheel was then installed in 
accordance with the manual. Extreme care was taken to 
select positions in the filter wheel which have no 
electronic filter compensation in the system. The manual 
denotes which positions these are. Once the system 
modification was complete, new calibration constants were 
computed and installed in the CATS system. Calibration 
constants are listed in Appendix D. 
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VII. THE NPS-NACIT SPLIT-FIELD POLARIMETER 



The split-field system approach possesses many 
advantages over external and internal filter placement 
normal to the path of the radiation, particularly in 
elimination of time lags and registration problems between 
orthogonal polarizations. The Naval Academic Center for 
Infrared Technology (NACIT) at the Naval Postgraduate School 
(NPS) has constructed and tested a split-field front end 
optical system for input to thermal imaging systems, 
specifically matched to the AGA-780 and the AMBER 
Engineering-4128. The Split-Field Polarimeter is the result 
of a conceptual design by Professor Alfred W. Cooper. The 
actual construction and design optimization was performed by 
Professor Eugene C. Crittenden and LT Peter H. Heisey, USN, 
as assistant. The general features of the NPS-NACIT Split- 
Field Polarimeter are: 

1. Simultaneous display and recording of orthogonally 

polarized images. 

2. Utilization of entire field-of-view for display. 

3. Adaptability to different imaging systems. 

4. Dual-band (8-12|im and 3-5nm) operation using 
interchangeable filters. 

5. F/1.87, 99mm focal length optics. 
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As mentioned in Chapt. IV for other systems with split-field 
optics, the 'Narcissus spot' is avoided by having the filter 
tilted relative to the path of the radiation. A background 
plate is also needed to give a uniform reflected radiation 
source into the images. The specifics of the system and its 
components and design will be discussed. 

A. SYSTEM DESIGN 

The split-field optical system consists of input and 
output compound lenses, polarizing splitter, and front- 
surface aperture limiting mirrors. The components are 
arranged in a vertical plane with respect to the ground. 

The system takes the collimated input radiation, focuses it 
down as it is split into two orthogonally polarized 
components, and folds and recombines the two split halves of 
the image while recollimating the light for processing by 
the thermal imaging system. In this process, the optical 
aperture is redefined so that only the scan lines of the 
central half of of each image are retained. The remaining 
3.5 x 7 degree halves are combined in adjacent halves of an 
intermediate 7x7 degree image plane. The 7x7 degree 
field of view of the optical system is thus split into two 
channels with orthogonal polarizations. Each channel is 
aperture limited into a 7x3.5 degree field to compose half 
of the final image. The original radiant flux in the image 
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is preserved in this way. An overall picture of the system 
in use with the AGA-780 Thermal Imaging System is shown in 
Fig. 11. In this setup, the entire imaging system is 
mounted on one lightweight mounting plate. Fig. 12 shows 
the inside of the box with all components exposed. Aluminum 
was used for all mountings and plates because of low weight, 
corrosion resistance, and easy machining. 




Figure 11. Top View of AGA with Split-System. 
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Figure 12. Side View of AGA with Split-System. 




Figure 13. Front View of Split-Field input 
lenses with AGA in background. 
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availability, corrosion resistance, and easy machining. The 
front view in Fig. 13 shows the exposed input lens. The 
long wave input lens for the AGA system is hidden by the 
backside of the system. The short wave lens not being used 
by the AGA is to the side and may be slid into place behind 
the optics when a short-wave polarizer is installed. Figure 
14 shows the image orientation and focal plane within the 
split-field optics. 




The incoming rays encounter two Germanium lenses (labeled 
lens system) and are refracted to a focal plane. In the 
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process, the polarizing splitter splits the beam into two 
paths with one beam vertically polarized with respect to the 
ground (transmitted beam) and the other reflected with 
horizontal polarization. The transmitted beam is reflected 
to a focal plane at the edge of a small mirror used by the 
reflected beam to be turned in the same direction (Fig. 14) . 
At this point, both beams have been recombined in a single 
focal plane. The beams then diverge and are reflected by the 
last mirror into two lenses identical to the first lens 
system (labeled lens system) . The radiation emerges 
afocally so that the AGA imager, positioned immediately 
behind the system can process the image for display. 

B. SPLIT-FIELD OPTICAL COMPONENTS 

The split-field system consists of 3 major components 
as stated in the previous chapter. The specific 
characteristics and operation of these components will now 
be addressed. 

1 . Lenses 

The lenses used by the system to generate the split 
image for input to the thermal imaging system consist of two 
identical assemblies each containing two germanium anti- 
reflection coated lenses. One assembly is used for the input 
and the other for the output side. The output assembly is 
reversed, end for end, relative to the input assembly. 
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The inner surface of the first lens of the input assembly 
and the inner surface of the last lens of the output 
assembly are aspheric. This reduces aberrations and provides 
a minimum spot size of 1.1 mrad, within the resolution of 
the AGA imager. The lens surfaces are coated for 

antireflection in both the 4|im and 12(om wave bands. This is 
necessary because the high index of refraction of germanium 
(4.0) would produce large reflection losses and scattered 
radiation if the surfaces were left uncoated. The overall 
focal length of each of the two lens assemblies is 99.2mm 
with an f number of 1.87 

2. Polarizing Splitter 

The most important component of the split system from 
the polarizing standpoint is the Al/Ge polarizing splitter. 
The custom designed splitter consists of a germanium 
substrate with an aluminum grid deposited on one side. The 
grid consists of .25(rm wide aluminum strips separated by 
.40nm for an overall period of . 65|im. The substrate is a 
round disc, 50mm in diameter x 2mm thick. Both sides of the 
splitter areanti-ref lection coated for the 8-12|im band. It 
is necessary to replace the splitter for the short-wave IR 
band, since the vendor found it impossible to coat the 
element for both bands. The same analysis done for the KRS- 
5 polarizers in Chapt. V was performed for the Ge splitter. 



49 



Data provided by the manufacturer (Graseby-Specac) is valid 
for normal incidence of the optical beam on the splitter. 
Values for K lf K 2 , and P are .81, 0.0, and 1.0, respectively 

at 10pm. This means the splitter passes 81% of the linearly 
polarized light perpendicular to the grid orientation and is 
~100% polarized. The manufacturers performance curves are 
provided in Appendix A. Extensive analysis was done to 
verify operation of the splitter; this will be addressed in 
section C. Results showed that there is a wavelength 
dependency of the radiation fraction that is reflected and 

transmitted. It is found that wavelengths between 7pm and 
9pm favor reflection and the opposite is true above 9pm. 

This is due to inefficiencies in the polarizing grid and 
anti-reflection coating. Results obtained by system 
calibration show this and will be discussed later. 

3. Gold Mirrors 

The folding operation of the optical system is 
performed by four gold first-surface ’protected' mirrors. 

The mirrors were purchased from Edmund Scientific Co. and 

are advertised to be 1/10 X flat at .5pm. Gold was chosen 
because of its high reflectivity, and thus, low emissivity. 
Figure 15 shows a plot of the emissivities at 9.92 pm for 
gold at a range of incidence angles. 
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Figure 15. Emissivity of Gold from eq. (9) . 



As expected, significant polarization effects occur only at 
angles close to grazing incidence. Most reflecting angles 
of the mirrors are near 45 degrees, so these effects are 
small (<1%) , but must be considered. 

C. SYSTEM CHARACTERISTICS 

The combined effects of the presence of the components 
in the optical path and mechanisms by which they perform 
their functions must be quantified to predict overall system 
performance. Irregularities in refracting and reflecting 
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surfaces, polarized reflection and emission, and system 
background radiance are effects which influence the measured 
radiance and resolution of the target. Experiments have 
been performed to characterize these effects to assist in 
predicting final image quality. 

1 . Resolution 

To obtain the full resolution possible by the AGA-780 
when used with the split-field system, effects induced by 
imperfections in the optics of the split-system must be 
minimized. Upon delivery of the splitter, visual resolution 
was tested using a standard bar chart. Reflection of the 
bar chart off the surfaces of the splitter showed no loss in 
resolution from direct viewing of the chart. It follows 
that imperfections in the surfaces would induce less effect 
in the longer wavelengths and would not be considered. 

A test was also done on the combined AGA-split system 
to discern if resolution was degraded by the presence of the 
split-field system. A sharp edged (cold) plate was placed 
in front of a black-body thermal source. Using the AGA 
optics alone, measurements were taken of the number of 
pixels in the image that were required to resolve the sharp 
edge between hot and cold surfaces. Along the direction of 
image raster scan, the resolution is about 2-3 pixels. 
Perpemdicular to the raster scan the resolution is slightly 
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better. Parallel to the raster scan the resolution is 
apparently limited by the frequency band pass of the 
amplifier electronics and the response time of the detector. 
In the direction perpendicular to the raster scan, the 
resolution is slightly higher because the response times of 
the detector and amplifier are not involved. In both 
parallel and perpendicular scan directions, no loss of 
resolution occurred in the image when the split-field optics 
were added. 

2 . MTF 

Modulation transfer function tests were performed for 
the AGA system with and without the split-field. The cutoff 
frequency for the AGA in the direction of raster scan is .3- 
.5 cycles/mrad which gives a resolution of 2-3mrad. The AGA 
system manual gives a system resolution of l.lmrad. This is 
apparently the geometrical resolution of the AGA optical 
system. The MTF data show that the overall system has a 2-3 
mrad resolution. There was no loss in spatial resolution 
with the split-field system; 

3. System Polarization Performance 

The main driver in the performance of the system is the 
polarizing splitter. Establishing numerical values for 
transmissivity, reflectivity, and absorption of the splitter 
for polarized and unpolarized target radiation requires 
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measurement of all radiance sources throughout the image and 
split-field system. 

a. Spectral Polarizer Performance 
The manufacturer of the polarizer provided 
spectral transmittance data for the polarizer at normal 
incidence to a reference beam. For use in the split system, 
it is necessary to know the spectral transmittance and 
reflection for 45 degree incidence. 

To measure the transmittance, the polarizer was placed 
at 45 degrees to the incoming beam in a Perkin-Elmer Model 
337 Infrared Grating Spectrophotometer. The photometer 
measures the characteristics in the same manner as the 
spectrophotometer used by the manufacturer. Because of 
polarization induced in the beam by the grating, the method 
described by Young et al . (Ref 13) , is used to compute the 
transmittance. This is the approach outlined in Chapt . VI. 
Appendix B shows the plots generated by the polarizer at 
normal and 45 degree incidence. Agreement with the 
manufacturers values are within 3% for normal incidence. 

For wavelengths less than 10.5fim, the transmissivity is 
higher than at normal incidence. This is probably due to 
the fact that the effective grid spacing is reduced by 
cos (45), when the grid is perpendicular to the plane of 
incidence of the beam, thus moving the curve left to shorter 
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wavelengths. It is necessary to mention that placing the 
filter in such a manner as to reflect ambient radiation into 
the detector does not contaminate the readings. The 
intensity radiated by a typical room temperature object 
compared to that of the IR source in the device is on the 
order of 0.001. 

Obtaining performance for reflection polarization 
required altering the conventional configuration of the 
polarizer to the beam to utilize the front surface as the 
effective transmitter. Figure 16 shows the arrangement used 
to obtain this. 




Figure 16. Configuration for polarization 
reflection measurement using the Perkin-Elmer 
M337 Infrared Grating Spectrophotometer. 
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Data for the linear polarization reflection efficiency is 
provided in Appendix B. Reference 13 provides a derivation 
of the measurement technique to obtain the transmittance of 
a polarizing grid. The difference in finding the efficiency 
here is that reflection is substituted for transmission in 
the equations. The radiant intensity with the polarizer 
removed was obtained by placing an aluminum mirror in place 
of the polarizer and resetting the 100% transmission 
indicator. The performance plot in appendix B shows this 
trace across the top of the graph. The designation for the 
reflected efficiency is K 1R . The equations are analogous to 
those described in chapter VI. It was not possible to 
measure K 2R (orthogonal polarization reflection) because of 
the unique arrangement and nonavailability of another 
similar polarizer. 

The degree of reflected polarization perpendicular 
to the grid, can be described as follows. Any radiation 
impinging on the grid with E vector perpendicular that is 
not transmitted must be reflected due to fact that the grid 
width is small compared to the wavelength, precluding 
current generation. This reflected component is the 
complement to the transmitted fraction of the radiation. 

This theory will be demonstrated in the following section. 
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b. Polarizer Performance Utilizing Radiance 

Analysis 

(1) Theoretical Representation 

An alternative method for measuring the fraction 
of polarized radiation transmitted and reflected by the 
polarizing splitter is to take radiance readings of the 
radiation passed with a known source polarization 
orientation. However, complications arise with this method 
when target temperature is near the temperature of the 
optical system, because the background irradiance that 
enters the detector is of the order of the target 
irradiance. To discriminate the target irradiance from the 
background, a mathematical description of the background 
effect is required. 

For the case of unpolarized radiation incident on 
the optical system, the following equations apply: 



I 

I 



T UUPOL 7, | T ± -L 

„ =f B I+f B I+S 

*ONPOL *| I R ± X 



( 19 ) 



where 



€j, r e R —the respective background emittances 



measured for each channel, independent 
of target irradiance 
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J^and J| represent the components of target irradiance 

perpendicular and parallel to plane of incidence of the 
polarizing splitter. These are what is actually measured by 
the thermal imaging system. For the entire analysis the 
convention used is that the polarizing splitter grid lines 
are perpendicular to the plane of incidence. The s T b,rb terms 
represent the radiant exitance of all the components of the 
system when the incident intensity into the system ~0. This 
is simulated by directing the optics at a pool of liquid 
nitrogen. The subscripts T and R designate the transmitted 
and reflected channel through the split-field optical 
system. T denotes the radiance fraction which passes 
through the polarizing splitter and is displayed as the top 
half of the split image; R denotes the other channel 
reflected by the splitter. The perpendicular and parallel 
subscripts describe the orientation of the E vector for the 
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irradiance relative to the plane of incidence of the 
splitter. For example, unpolarized radiation is commonly- 
represented by 



J = ( J /2) + ( J /2) (20) 

To eliminate unknown terms in Eq. (19), measuring 
a purely linearly polarized source in two orthogonal 
orientations yields a simpler set of equations. For the 
first case of purely polarized target irradiance 
perpendicular to the plane of incidence (perpendicular to 
the grid) , 




I,+S„ +2, 




J. + 2„ +2 , 



( 21 ) 



The f TjR terms denote that fraction which goes into the 
respective channels from splitter reflection and 
transmission performance. The target consisted of a 
polarizing filter in front of a black body source; 
therefore, there is some filter emittance that must be 

considered. The added e F terms represent the emittance of 
the polarizing source for each channel, for each orientation 
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of the target producing filter. It can be seen that there is 
a fractional contribution to each channel to be measured T 
and R. This is because it is known that the transmittance 
of the grid is not 1.0. For the case where the target 
radiance is plane polarized perpendicular to the plane of 
incidence (parallel to the grid) , 



J 






J +E 




( 22 ) 



There are few terms seen in the T channel because it is 
known from the spectral performance curves that the 
transmitted E components are ~100% polarized. Equation 21 
shows the intensities seen in the two channels with a 
vertically, parallel (p) linearly polarized source and eq. 

22 with 90 degree rotation of the target E vector, s 
polarized. 

(2) Experimental Radiance Measurement 
To compute the performance of the system using the 
equations in the previous section, it is necessary to take 
careful radiance measurements of each term in the equations 
with a linearly polarized target. Appendix C contains the 
values obtained in the laboratory for these measurements. 
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To compare the performance of the polarizing splitter 
channels, the values in Appendix C are plugged into eqns . 
(21) and (22) and then subtract the equations from each 
other. The relationships obtained are 

(f-f)I= 21.9 

T \ 1 

f D I = 22.4 

A. J - 



Measurements showed that there is no dependency of target 
intensity on polarization orientation, so T ( = . 

Therefore, 



<VV « no 

!_= 0. 98 



(23) 



This equation allows calculation of the reflected fraction 
of the incident radiation which is normally passed through 
by the grid. The implications of Eq. (23) are that the 
transmittance of linearly polarized irradiance perpendicular 
to the grid is always greater than the reflection efficiency 
of the orthogonal component from the grid front. This is 
seen with the performance curves in appendix B. It also 
allows correlation of splitter performance using the 
radiance method with that seen in the spectral performance 
curves. To compare the radiance method, it will be 
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assumed that the non-transmitted E vector perpendicular to 
the grid undergoes Snell reflection and f T +f R =l. The 
spectrophotometer data measured in the previous section 

shows that at 9.52[im, the above expression = 0.97. 

Agreement is seen between the two methods at that 
wavelength, but the ratio varies strongly at other 
wavelengths. Schematic representation of what happens at 
the splitter within the optical system is shown in Figs. 17- 
19. Figure 17 indicates the fractional performance in the 
LWIR for unpolarized input. For this case, 58% of the 
background is reflected into the transmitted (T) channel and 



Reflected conpent 
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Background Plate 



Figure 17. Schematic of polarizing splitter 
for unpolarized incident radiation. Grid 
perpendicular to plane of incidence. 
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combined with 41% of the target irradiance. The reflected 
(R) channel combines 42% of the background with 59% of the 
target irradiance. Figure 18 indicates the fractional 



R*£l«Cw«d Coopon«ot 




1 

Background Plate 



Figure 18. Schematic of polarizing splitter 
for horizontally polarized incident 
radiation. Grid perpendicular to plane of 
incidence . 



performance in the LWIR for horizontally polarized input. 
For this case the transmitted beam contains 58% of the 
irradiance from the background plate and essentially no 
target radiation. Figure 19 shows the equivalent situation 
for vertically polarized input, which shows 83% of the 
target irradiance with 58% of the background irradiance in 
the transmitted (T) beam, and 17% of the target irradiance 
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Reflected component 




Figure 19. Schematic of polarizing splitter 
for vertically polarized incident radiation. 
Grid perpendicular to plane of incidence. 



with 42% of the background plate irradiance in the reflected 
(R) beam. 

The radiance analysis of the polarizing splitter 
performance provides a gross analysis and verification of 
the data provided by the manufacturer and by the grating 
spectrophotometer. This representation of transmittted and 
reflected target irradiance must be considered when 
processing split-field images. 
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VIII. IMAGE ANALYSIS 



To prove the concepts presented and to demonstrate the 
effects of polarized filtering on image contrast, it is 
necessary to perform a preliminary image analysis on images 
recorded with the thermal imaging system in two 
configurations; internal polarizing and split-image optical 
system polarizing. Two separate experiments were conducted 
to demonstrate system performance. 

In March 1995, 1338 images were taken during a field 
experiment with the R/V POINT SUR from the Monterey Bay 
Aquarium Research Institute marine operations building at 
Moss Landing, CA. All of the polarization filtering was 
done internally. In August of the same year, images were 
recorded using the NACIT Split-field Polarimeter from the 
Naval Postgraduate School beach property of the USS John 
McCain, an Arleigh Burke Class DDG. The images recorded 
show contrast improvement effects from polarization 
filtering. 

A. DEFINITION OF CONTRAST AND TARGET -BACKGROUND RADIANCE 

DIFFERENCE 

When an observer views an image, the critical parameter 
for detecting and recognizing the target amidst the 
background is target contrast. Numerous definitions are 
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currently used to define contrast. Common examples of these 
are [23] 



AT 



target background 

AT “ 

background 



N 



Contrast- 



Contrast= targec back ? round 



(24) 



N 



target background 



Here N represents radiance. It is desirable to compare the 
contrast of images which may be separated in time. In this 
case, changes in system gain cause the apparent digitized 
image values to change for similar viewing conditions. 

Also, changing instrument settings changes N for every image 
when the controls are altered. Therefore, when target and 
background conditions vary little between images, the 
difference in N may be a more accurate indicator of 
contrast. Differences in filter transmittance cause the 
conventional contrast equations to be in error so that any 
artificial thermal level changes induced by this are ignored 
for filters with the same calibration curve slope. This is 
the case for the internal filters (Appendix E) . To 
accommodate this, the analysis will only consider 
differences in the apparent radiances of the target and 
background. The contrast parameter to be measured is termed 
Target-Background Radiance Difference (TBRD) . This is the 
numerator of both contrast equations in Eq. (24) . 
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TBRD = A N = <N. 



(25) 



„ > “ <N, . > 

target background 



The average values of N are used because of the varying 
composition of the image with target distance, system noise, 
and measurement points in the image. If the image radiation 
is polarization filtered, the resultant TBRD is a function 
of the amount of filtering only if all optical components 
and scene characteristics remain constant. It has been 
shown that many different components of an image can display 
different degrees of polarization. If background irradiance 
is suitably filtered more than target irradiance, TBRD will 
increase . 

When comparing two images, it is useful to quantify the 
increase in image contrast improvement. The TBRD 
improvement factor is defined by 



AN 

TBRD. =AN. = 

improvement imp 



(26) 



B. FIELD MEASUREMENTS WITH INTERNAL FILTERS 

From May 16-23, 1995, polarized thermal imaging 
experiments were conducted from the Monterey Bay Aquarium 
Research Institute using the R/V POINT SUR as target. 
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1338 images were recorded using internally filtered and 
unfiltered images in the 8-12|rm wave band. The ship 

presented all aspects with ranges varying from .25 nmi to 
the range limit of the day. Appendix E carries a complete 
listing of the images and a small picture of each with the 
filename shown on the bottom of each image. Vertically and 
horizontally polarized images will be compared for target- 
background radiance difference. 

1 . Visual Image Comparison 

Selected pairs of images demonstrate the visual TBRD 
and contrast improvement between vertically and horizontally 
polarized images. Unpolarized images were not considered 
for this comparison because of the nature of the data 
analysis process, and the change in system settings 
necessary to obtain useable images due to the changing of 
system gain and altering the byte level of the images 
significantly. The definition of TBRD would not be valid 
for comparison to the polarized images and further 
compensation to TBRD would be necessary. The images are 
selected from the compilation listing in Appendix E. The 
contrast between the two orthogonally polarized members of 
each pair of images in Fig. 20a and 20b shows the lower 
apparent radiance of the sea water background in the 
horizontal polarization. It is necessary to mention that 
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the images were printed on a Hewlett-Packard LaserJet 4M 
600dpi printer. Differences in the images are more easily 
seen in a grayscale color table on an SVGA monitor. The 
bright area seen in the top right corner of each image is a 
constant temperature source created by inserting a mirror in 
the field of view, reflecting a calibrated black body 
source. This provides an unpolarized reference in each 
image to compare filter transmittance and provide 
temperature calibration confirmation. The source is not 
discussed in this analysis. 
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M5L16051 16May95 
Polarization: Horizontal 
TBRD=108 . 5 




M5L17195 17May95 
Polarization: Horizontal 
TBRD=58 . 3 




M5L16052 16May95 
Polarization: Vertical 
TBRD=73 . 2 
TBRD^l.4 8 




M5L17196 17May95 
Polarization: Vertical 
TBRD=44 . 0 
TBRD imp =l . 33 



Figure 20a. Comparison of vertically and horizontally 
polarized LWIR Band (8-12jam) images of RV POINT SUR. 



70 




M5L23014 23May95 
Polarization: Horizontal 
TBRD=22 . 8 




M5L17136 17May95 

Polarization: Horizontal 

TBRD=38 . 4 




M5L23015 23May95 
Polarization: Vertical 
TBRD=10 . 2 
TBRDimp=2 .24 




M5L17137 17May95 

Polarization: Vertical 
TBRD=32 . 4 
TBRD^l.18 



Figure 20b. Comparison of vertically and horizontally 
polarized LWIR Band (8-12pm) images of R/V POINT SUR 
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2 . 



Numerical TBRD Measurement. 



The images presented in Fig. 20 were extracted groups 
of images with the same aspects, ranges, and path 
conditions. The average radiances are calculated by 
averaging over the target and background areas within an 
area around the ship. Appendix F contains the IDL code to 
perform this operation. To compare the TBRD for a number of 
images, the images must be taken with target and background 
conditions as constant as possible for trends in the data to 
appear. Figures 21-23 show TBRD measurements made on 
different sets of images grouped by ship range and aspect. 
The vertical axis gives the TBRD in dimensionless units. 

The 8 bit A/D converter sets the maximum thermal level 
reading based on system gain to 255. This irradiance 
reading is independent of system calibration settings and 
measures the raw data from the detector with no compensation 
for temperature accuracy. The horizontal axis gives range 
computed from ship image height measured in pixels. Smaller 
targets appear to the right side of the plot, which shows 
TBRD as a function of range. The general trend of TBRD 
reduction with target size is due to increased path 
attenuation and path radiance with range. The legend on the 
plots denotes the polarization and whether TBRD is 
calculated above or below the horizon. The average TBRD 
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ML5230 10-030 TBRD vs. Range 
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Figure 21. TBRD Improvement-Starboard Aspect 
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Ship Range, Nautical Miles, from Ship Size 



M5L17180-200 TBRD vs. Range 
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Figure 22. TBRD Improvement-Port and Stb Aspect 
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Ship Range, Nautical Miles, from Ship Size 



ML516050-16070 TBRD vs. Range 
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Figure 23. TBRD Improvement-Bow Aspect 
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Ship Range, Nautical Miles, from Ship Size 



improvemement on each plot is calculated only for the below 
horizon case against the sea water background. It is 
expected that TBRD against the sky should be nearly the same 
in both polarizations. This data does not contain any 
meteorological, target temperature, or detailed ship 
position data which is available for analysis improvement. 

Very general trends are seen for TBRD improvement, but 
are not as apparent as the average improvement calculation 
numbers in the plot. Ship size varies greatly in the plots 
due to variation in the numbers of pixels recorded in the 
vertical direction and is most pronounced in the images 
taken in the near ranges. It is seen that there is some 
polarization effect in the above horizon TBRD. This is 
probably due to horizon selecton in the image and possible 
polarization of atmospheric radiance close to the horizon. 
Analysis was not performed on the polarization of the sky. 
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C. FIELD MEASUREMENTS WITH SPLIT-FIELD POLARIMETER 

Preliminary test and evaluation imaging was done on 3 
August, 1995 with the split-field polarimeter optical 
system. The target is the USS John McCain, an Arleigh Burke 
class destroyer present in Monterey Bay for a port visit at 
the Naval Postgraduate School. Thirteen split-field and 10 
test images were taken of the ship and black-body calibra- 
tion sources. The target was at a stationary range of 
1.85km, determined by target geometrical size. All images 
are port aspect. All of the images taken are listed by 
filename in Appendix G. Because of the similarity in image 
composition, only one image will be analyzed in depth. 

Figure 24 shows the image chosen which is optimized for 
focus, level gain setting, and system adjustment. Images 
are recorded inverted by the AGA and this image has been 
artificially flipped upright for viewing. The top image is 
the reflected component of the polarizing splitter, 
horizontally (s) polarized, and the bottom is the trans- 
mitted vertically (p) polarized image. Consistent with 
results seen for the internally polarized images is the 
increase in contrast of the top image. Discussion of 
specific composition of each channel is in Chapter VII. 
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Figure 25 shows a vertical profile taken of the image 
through the hot stack. From the known geometry of the ship, 
no significant polarization is expected for the stack which 
is just aft of the large notch in the top of the 
superstructure. Polarization of emission could exist in the 
ship panels on the bow as they slant outward toward the 
horizontal with height. The bottom of the profile 
corresponds to the bottom of Fig. 24. The x-axis denotes 




Figure 24. NACIT Split-field polarized image, 
USS John McCain, 3Aug95 
Top-Vertically Polarized TBRD=168 

Bottom-Horizontally Polarized TBRD=200 
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the 8 bit digitized value of radiant intensity. TBRD for 
the top image is 200 and 168 for the bottom. This gives a 
TBRD improvement factor = 1.19. This is with reference to 
the average value of the sea water background taken in the 
foreground along the azimuth of the stack in the middle. 

The fog above the image is assumed unpolarized. There is an 
overall increase in image intensity in the vertically 
polarized image due to the behavior of the splitter and the 
background plate (Ch 7) . 
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Figure 25. Vertical Profile of Split Image, M0803011 
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The added reflectance from the back of the splitter raises 
the overall level of the vertically polarized image. This 
is another justification for using TBRD to compare image 
contrast differences. 

1 . Digital Image Subtraction 

The goal of the split system is to provide simultaneous 
images of the two polarizations which allows digital 
subtraction of the images to visually represent TBRD 
improvement. An IDL program was written to accomplish a 
pixel by pixel overlay of the images of the ships (Appendix 
F) . The program finds the 'hottest' pixel in the image and 
aligns the pixel addresses in the top and bottom while 
subtracting the values. The resulting image displays a 
value of 0 (black) where the images had the same intensity. 
The areas in the image where polarization effects occur are 
the brightest. Figure 26 is an image of the top channel 
minus the bottom channel compensated for background 
contribution to the bottom image. This was done by 
adjusting the threshold of intensity for the bottom image. 
This essentially removes any polarization effects seen on 
the ship because of their small effects. 

It is clearly seen that the water shows the largest 
amount of difsference in the image. Irregularities in 
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Figure 26. Digitally Subtracted top and bottom 
split image channels. M0803011, 3Aug95. 

in the image plane show on both sides of the subtracted 
image. Techniques are under development to balance the 
intensity of the two images, remove background 
irregularites, and improve alignment of the two images. 
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but these techniques were not available at the time 
these images were processed. The images are shown here to 
illustrate the general utility and feasibility. 
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IX. DISCUSSION AND CONCLUSIONS 



The split-field polarizing system designed and con- 
structed in this work, demonstrates the utility of simulta- 
neous image storing for processing, while preserving image 
resolution. 

Subjective and numerical analysis of the data from 
images taken with the split-field system demonstrated excel- 
lent performance in image quality, contrast improvement, and 
utility. Polarization filtering performance for a sea water 
background, for both the split-field system and for direct 
viewing with interchangeable filters in the AGEMA imager, 
gave results comparable to those obtained in the MAPTIP 
experiments where an external BaFl substrate polarizer was 
used. This was determined by comparing contrast calculated 
in terms of TBRD from previous measurements with values from 
the split-field and internal filtering systems. Both the 
split-field and internally filtered direct viewing systems 
eliminated the narcissus effects that degrade the use of 
external filters. 

Observations of the target scenes show quantitatively 

that : 

1. Digital subtraction of simultaneously recorded split- 
field images is a viable technique for contrast improvement. 
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2. In the 8-12 |im region, vertical polarization emission of 
the sea surface is the predominant component for filtering. 

3. Split-field polarized filtering gives preliminary TBRD 
improvement of contrast of 1.19. This is driven by the 
splitter performance and background plate effects. 

4. Direct viewing with internal polarization filtering is 
capable of 1.2- 1.5 TBRD improvement of contrast for sea- 
water emission, but digital subtraction is not possible 
because of time delay between acquisition of the 
orthogonally polarized images. 

5. Both the split-field and direct viewing internal filter- 
ing techniques provide contrast improvement on the order of 
that seen in the MAPTIP experiments. 

6. The combined system designed and constructed for this 
project provides dual use of the AGA-780 with the split- 
field system and as a direct-view internally-filtered sys- 
tem. This provides two highly effective polarizing systems 
for use in one location with flexibility in system selec- 
tion . 
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X. RECOMMENDATIONS FOR FUTURE WORK 



Further work needs to be performed in all aspects of 
polarized imaging using the split system and internal 
polarizing systems. In March, 1996 an Electro-Optical 
Propagation Assessment in the Coastal Environment (EOPACE) 
experiment phase will take place near San Diego CA. In- 
depth demonstration and field implementation of both 
techniques will take place at that time. Specific 
concentration areas for both systems in the intermediate 
time frame include: 

Internal Polarizers: 

-Quantify the spectral transmissivity experimentally 
for each polarizing filter to verify percent polarization 
performance . 

-Re-calibrate the AGA system with the filters installed 
prior to use to give level equalization for both filters. 

Split-image Polarimeter: 

-Conduct laboratory experiments to explore and remedy 
system induced gradient across the field-of-view. 

-Realign Split-field optics to eliminate any 
vignetting, focussing differences between the two channels, 
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and edge effects due to optical surface misalignment. 

-Install a temperature control system on the background 
plate to maintain constant radiance contribution to the 
image for varying environmental conditions. 

-Insert a small temperature source in the intermediate 
focal plane of the split-field optics to give a constant 
calibration source for real time split channel comparison. 

Overall AGA-780 Thermal Imaging System: 

-Install the new NEC VERSA P/75 75MHz Pentium Laptop 
computer as the display and file storage server to replace 
the current CATS 386 computer. 

-Install the new 12 bit A/D converter for increasing 
system dynamic range and allowing higher data resolution for 
analysis . 
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APPENDIX A. 



PERFORMANCE CURVES FOR THE POLARIZING FILTERS 



This appendix contains the performance curves for the 
polarizing filters used in the internal and split image 
configuration. The interpretation of the data is included in the 
text of the appropriate chapters. 

The data were provided by the manufacturer of the 
polarizers, Graseby-Specac . The measurements were taken with a 
spectrometer using the actual filter shipped. Where two filters 
were required for the crossed-grid measurement (trace 3) , another 
polarizer nearly identical to the one shipped is used. 

Trace 1- transmittance with grid perpendicular to the E vector of 
incident radiation. 

Trace 2- transmittance with grid parallel to the E vector of 
incident radiation. 

Trace 3- transmittance of two filters with grids crossed. 

Vertical Axis- wavenumber/wavelength 
Horizontal Axis- transmittance 
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Figure A2 . Performance curves for KRS-5 9.5mm polarizer (#75) 
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APPENDIX B. 

SPECTRAL TRANSMITTANCE /REFLECTANCE OF THE 50 MM POLARIZER 

This appendix contains the Perkin-Elmer 
spectophotometer traces for the 50mm polarizing splitter. 

In Fig. Bl, the traces with the x designation are for the 
polarizer oriented 45 degrees to the incoming beam. The K 1N 
trace is the combination of the traces and gives the overall 
transmittance of the splitter normal to the incoming beam. 
The x symbols along the K 1N trace are the overall 
transmittance of the splitter at 45 degrees to the incoming 
beam. It is seen that the peak occurs at a slightly shorter 
wavelength than at normal incidence. 

In Fig.B2, K 1R is the overall reflection efficiency in 
the configuration described in Chapt . VI. 
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APPENDIX C. 



MEASURED RADIANCE VALUES FOR THE NPS-NACIT SPLIT-FIELD 

POLARIME TER 



The following data were taken as input to the radiance 
analysis of polarizing splitter performance. 

Lab Conditions: 

Room Temp 17C, 

Source-Agema Black-Body 660C 

Polarized Source-BaFl IR Polarizer Pol%~99% 

Sensor-AGA 780 LW Band, f/1.87, Range 2 
Background Plate Temp- ~0C (heavy frost surface) 
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APPENDIX D. 



CALIBRATION CURVES FOR FILTERED AND UNFILTERED AGA-780 
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LW,No filter (position 6) 




ONtavaa na/an tvwh3Kl 



Figure D-l. Calibration for the AGA with no filters 
present . 
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ONIOVBy 13A3H TVWH3H1 



Figure D2 . Calibration Data for AGA with internal filter. 
Pos 5 denotes one of the small internal filters located in 
position 5 of the AGA filter wheel. Position 4 is the other 
with orthogonal filter orientation. A small difference is 
seen in thermal level with temperature between the two 
filters . 
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0NI0V3M T3A31 TVWiGHJ. 



Figure D3 . Calibration curve for the Split-field 
Polarimeter. This shows level reading with temperature for 
the two channels of the split-field polarimeter. LWPOL is 
for the large BaFl external polarizer for comparison. Upper 
spot denotes the splitter transmitted channel imaged by the 
AGA. 
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APPENDIX E. 

COMPILATION OF MOSS LANDING IMAGES 
MAY, 1995 



The following is a complete image compilation of the entire 
Moss Landing experiment. The data tabulated is incomplete and 
only shows what is necessary for preliminary TBRD comparison of 
images. A sheet depicting a number of representitive images is 
appended to this list. Images corresponding to each of the table 
entries are available in digital form on hard disk at the NACIT 
laboratory. 
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The title header gives data about each of the images taken. 

Pol-polarization orientation of the image, h, v-horiz, vert . x-image not 
analyzable for one reason (see image) . u-unpolarized image. 

Level-thermal level setting of the AGA control knob. 

Range-thermal level range setting of the AGA control. 

Offset-Offset of the calibration constant for image temperature compensation. 
Ovrflw, Undrflw-number of pixels in the image with radiance values above or 
below the 0-255 range. 



Filename 


date 


time 


pol 


level 


range 


offset 


ovrflw 


undrflw 


M5L17001. IMG 


5/16/95 


12:43:06 


X 


44.70 


20. 


2. 


0 . 


0 . 


M5L17002 . IMG 


5/16/95 


12:43:44 


X 


44.70 


20. 


2. 


0 . 


0 . 


M5L17003 . IMG 


5/16/95 


12:45:09 


X 


44.70 


20. 


2. 


0 . 


0 . 


M5L17004 . IMG 


5/16/95 


12:45:20 


X 


44.70 


20. 


2. 


0 . 


0 . 


M5L17005. IMG 


5/17/95 


9:11:07 


h 


46.90 


2. 


0 . 


246. 


0 . 


M5L17006. IMG 


5/17/95 


9:11:20 


V 


46. 90 


2. 


0 . 


294. 


0 . 


M5L17007. IMG 


5/17/95 


9:11:59 


u 


44.40 


5. 


0 . 


202. 


0 . 


M5L17008 . IMG 


5/17/95 


9:12:20 


u 


46.40 


5. 


0 . 


4. 


16482. 


M5L17009 . IMG 


5/17/95 


9:13:09 


h 


46.60 


2. 


0 . 


297. 


0 . 


M5L17010 . IMG 


5/17/95 


9:13:20 


V 


46.60 


2. 


0 . 


410. 


0 . 


M5L17011. IMG 


5/17/95 


9: 13:51 


V 


47.40 


2. 


0 . 


210. 


0 . 


M5L17012 . IMG 


5/17/95 


9:14:15 


u 


44.40 


5. 


0 . 


280. 


0 . 


M5L17013 . IMG 


5/17/95 


9:15:07 


u 


46.60 


5. 


0 . 


3. 


18368. 


M5L17014 . IMG 


5/17/95 


9:17:06 


h 


46.90 


2. 


0 . 


252. 


0 . 


M5L17015. IMG 


5/17/95 


9:17:16 


V 


46. 90 


2. 


0 . 


391. 


0 . 


M5L17016. IMG 


5/17/95 


9:17:33 


u 


45.40 


5. 


0 . 


100. 


3916. 


M5L17017. IMG 


5/17/95 


9:17:52 


u 


47.40 


5. 


0 . 


2. 


18902. 


M5L17018 . IMG 


5/17/95 


9:21:06 


h 


46.90 


2. 


0 . 


210. 


0 . 


M5L17019 . IMG 


5/17/95 


9:21:16 


V 


46.90 


2. 


0 . 


305. 


0 . 


M5L17020. IMG 


5/17/95 


9:21:38 


u 


44.90 


5. 


0 . 


169. 


0 . 


M5L17 02 1 . IMG 


5/17/95 


9:21:57 


u 


47.40 


10. 


0 . 


0 . 


0 . 


M5L17022 . IMG 


5/17/95 


9:22:55 


h 


47.40 


2. 


0 . 


77. 


978. 


M5L17023 . IMG 


5/17/95 


9:23:04 


V 


47.40 


2. 


0 . 


178. 


0 . 


M5L17024. IMG 


5/17/95 


9:23:21 


u 


44.90 


5. 


0 . 


186. 


0 . 


M5L17 02 5 . IMG 


5/17/95 


9:23:39 


u 


48.10 


10. 


0 . 


0 . 


8703. 


M5L17026 . IMG 


5/17/95 


9:25:13 


h 


47.10 


2. 


0 . 


161. 


0 . 


M5L17 027 . IMG 


5/17/95 


9:25:23 


V 


47.10 


2. 


0 . 


264. 


0 . 


M5L17028 .IMG 


5/17/95 


9:25:43 


u 


45.20 


5. 


0 . 


117. 


1705. 


M5L17029 . IMG 


5/17/95 


9:25:59 


u 


48.60 


10. 


0 . 


0 . 


9494. 


M5L17030 . IMG 


5/17/95 


9:27:01 


h 


46. 60 


2. 


0 . 


256. 


0 . 


M5L17031 . IMG 


5/17/95 


9:27:09 


V 


46. 60 


2. 


0 . 


343. 


0 . 


M5L17032 . IMG 


5/17/95 


9:27:28 


u 


44.90 


5. 


0 . 


154. 


1 . 


M5L17033 . IMG 


5/17/95 


9:27:42 


u 


48.60 


10. 


0 . 


0 . 


16600. 


M5L17034 . IMG 


5/17/95 


9:34:30 


h 


47.10 


2. 


0 . 


192. 


0 . 


M5L17035. IMG 


5/17/95 


9:34:38 


V 


47.10 


2. 


0 . 


300. 


0 . 


M5L17036. IMG 


5/17/95 


9:34:53 


u 


44.90 


5. 


0 . 


178 . 


0 . 


M5L17037 . IMG 


5/17/95 


9:35:06 


u 


48.10 


10. 


0 . 


0 . 


7669. 


M5L17038 .IMG 


5/17/95 


9:35:39 


h 


47.10 


2. 


0 . 


175. 


0 . 


M5L17039 . IMG 


5/17/95 


9:35:44 


V 


47.10 


2. 


0 . 


276. 


0 . 


M5L17040 . IMG 


5/17/95 


9:36:05 


u 


45.40 


5. 


0 . 


60. 


5958. 


M5L17041 . IMG 


5/17/95 


9:36:21 


u 


47.90 


10. 


0 . 


0 . 


5703. 


M5L17042 . IMG 


5/17/95 


9:37:06 


h 


47.10 


2. 


0 . 


31. 


4495. 


M5L17043 . IMG 


5/17/95 


9:37:14 


V 


47.10 


2. 


0 . 


163. 


249. 


M5L17044 . IMG 


5/17/95 


9:37:28 


u 


45.20 


5. 


0 . 


51. 


6364. 


M5L17045. IMG 


5/17/95 


9:37:42 


u 


47.90 


10. 


0 . 


0 . 


7616. 


M5L17046. IMG 


5/17/95 


9:42:58 


h 


46.40 


2. 


0 . 


277. 


0 . 


M5L17047 . IMG 


5/17/95 


9:43:07 


V 


46.40 


2. 


0 . 


398. 


0 . 


M5L17048 . IMG 


5/17/95 


9:43:19 


u 


44.70 


5. 


0 . 


176. 


19. 


M5L17049. IMG 


5/17/95 


9:43:33 


u 


48.10 


10. 


0 . 


0 . 


7633. 


M5L17050 . IMG 


5/17/95 


9:44:16 


h 


46.40 


2. 


0 . 


305. 


0 . 
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M5L17051 . IMG 


5/17/95 


9:44:25 


V 


46.40 


2. 


0. 


435. 


0. 


M5L17052 . IMG 


5/17/95 


9:44:41 


u 


44.90 


5. 


0. 


157. 


322 . 


M5L17 053 . IMG 


5/17/95 


9:44:54 


u 


48.10 


10. 


0. 


0. 


9403. 


M5L17054 . IMG 


5/17/95 


9:45:28 


h 


46.90 


2. 


0. 


189. 


0. 


M5L17055. IMG 


5/17/95 


9:45:36 


V 


46.90 


2. 


0. 


304. 


0. 


M5L17056 . IMG 


5/17/95 


9:45:51 


u 


44.20 


5. 


0. 


254. 


0. 


M5L17057 . IMG 


5/17/95 


9:46:03 


u 


48.10 


10. 


0. 


0. 


8709. 


M5L17058 . IMG 


5/17/95 


9:47:34 


h 


46.90 


2. 


0. 


222. 


0. 


M5L17059. IMG 


5/17/95 


9:47:43 


V 


46.90 


2. 


0. 


343. 


0. 


M5L17060. IMG 


5/17/95 


9:48:00 


u 


44.70 


5. 


0. 


207. 


0. 


M5L17061 . IMG 


5/17/95 


9:48:16 


u 


48.10 


10. 


0. 


0. 


9469. 


M5L17062 .IMG 


5/17/95 


9:49:00 


h 


47.10 


2. 


0. 


73. 


1898. 


M5L17063 . IMG 


5/17/95 


9:49:08 


V 


47.10 


2. 


0. 


238. 


0. 


M5L17064 . IMG 


5/17/95 


9:49:21 


u 


44.70 


5. 


0. 


240. 


0. 


M5L17065 . IMG 


5/17/95 


9:49:37 


u 


48.40 


10. 


0 . 


0. 


8762. 


M5L17066 . IMG 


5/17/95 


9:53:32 


h 


46.90 


2. 


0. 


370. 


0. 


M5L17067 . IMG 


5/17/95 


9:53:41 


V 


46.90 


2. 


0. 


522. 


0. 


M5L17068 .IMG 


5/17/95 


9:53:57 


u 


44.90 


5. 


0. 


211. 


0. 


M5L17069 . IMG 


5/17/95 


9:54:11 


u 


48.60 


10. 


0. 


0. 


12991 


M5L1707 0 . IMG 


5/17/95 


9:54:35 


h 


47.40 


2. 


0. 


254. 


0. 


M5L17071 . IMG 


5/17/95 


9:54:44 


V 


47.40 


2. 


0. 


300. 


0. 


M5L17072 .IMG 


5/17/95 


9:55:02 


u 


44.90 


5. 


0. 


192. 


100. 


M5L1707 3 . IMG 


5/17/95 


9:55:15 


u 


48.60 


10. 


0. 


0. 


17985 


M5L17074 . IMG 


5/17/95 


9:57:35 


h 


47.40 


2. 


0. 


182. 


0. 


M5L17075. IMG 


5/17/95 


9:57:46 


V 


47.40 


2. 


0 . 


307. 


0. 


M5L1707 6 . IMG 


5/17/95 


9:58:01 


u 


45.20 


5. 


0. 


196. 


0. 


M5L17 077 . IMG 


5/17/95 


9:58:14 


u 


48.40 


10. 


0. 


0. 


10557 


M5L17078 . IMG 


5/17/95 


9:58:37 


h 


47.10 


2. 


0. 


444. 


0 . 


M5L17079. IMG 


5/17/95 


9:58:47 


V 


47.10 


2. 


0. 


655. 


0. 


M5L17080 . IMG 


5/17/95 


9:59:08 


u 


44.70 


5. 


0. 


191. 


122. 


M5L17081 . IMG 


5/17/95 


9:59:27 


u 


47.40 


10. 


0. 


0. 


5246. 


M5L17082 . IMG 


5/17/95 


10:07:29 


h 


47.40 


2. 


0. 


268. 


0. 


M5L17 083 . IMG 


5/17/95 


10:09:43 


h 


47.90 


2. 


0. 


356. 


0. 


M5L17084 . IMG 


5/17/95 


10:10:02 


V 


47.90 


2 . 


0 . 


630. 


0. 


M5L17 085 . IMG 


5/17/95 


10:10:19 


u 


44.40 


5. 


0. 


288. 


0. 


M5L17086 . IMG 


5/17/95 


10:10:35 


u 


48.10 


10. 


0. 


0. 


6029. 


M5L17087 . IMG 


5/17/95 


10:11:16 


h 


47.60 


2 . 


0. 


704. 


0. 


M5L17088 . IMG 


5/17/95 


10:11:35 


V 


48.10 


2. 


0 . 


517. 


0 . 


M5L17089. IMG 


5/17/95 


10:12:01 


u 


44.40 


5. 


0 . 


316. 


0 . 


M5L17090. IMG 


5/17/95 


10:12:18 


u 


47.90 


10. 


0 . 


0 . 


6477. 


M5L17091 . IMG 


5/17/95 


10:14:11 


h 


47.90 


2. 


0 . 


447. 


0 . 


M5L17092 . IMG 


5/17/95 


10:14:19 


V 


47.90 


2. 


0 . 


724. 


0 . 


M5L17093. IMG 


5/17/95 


10:14:47 


u 


44.70 


5. 


0 . 


290. 


0 . 


M5L17094. IMG 


5/17/95 


10:15:20 


u 


47.90 


10. 


0 . 


0 . 


6163. 


M5L17095. IMG 


5/17/95 


10:15:57 


h 


47.90 


2. 


0 . 


461. 


0 . 


M5L17096 . IMG 


5/17/95 


10:16:09 


V 


47.90 


2. 


0 . 


617. 


0 . 


M5L17097 . IMG 


5/17/95 


10:16:29 


u 


45.20 


5. 


0 . 


202. 


3449. 


M5L17098 . IMG 


5/17/95 


10:16:44 


u 


48.40 


10. 


0 . 


0 . 


8100. 


M5L17099 . IMG 


5/17/95 


10:21:54 


h 


48.60 


2 . 


0 . 


354. 


0 . 


M5L17100. IMG 


5/17/95 


10:22:04 


V 


48.60 


2. 


0 . 


618. 


0 . 


M5L17 101 . IMG 


5/17/95 


10:22:20 


u 


45.20 


5. 


0 . 


414. 


0 . 


M5L17102 . IMG 


5/17/95 


10:22:37 


u 


48.10 


10. 


0 . 


0 . 


6539. 


M5L17103 . IMG 


5/17/95 


10:22:54 


h 


48.60 


2. 


0 . 


629. 


0 . 


M5L17104 . IMG 


5/17/95 


10:23:02 


V 


48.60 


2 . 


0 . 


885. 


0 . 


M5L17105 . IMG 


5/17/95 


10:23: 17 


u 


44.70 


5. 


0 . 


504. 


0 . 


M5L17106. IMG 


5/17/95 


10:23:35 


u 


48.80 


10. 


0 . 


0 . 


17517 


M5L17107 .IMG 


5/17/95 


10:26:59 


h 


48.60 


2. 


0 . 


991. 


0 . 


M5L17108 . IMG 


5/17/95 


10:27:22 


V 


48.60 


2. 


0 . 


913. 


0 . 


M5L17109 . IMG 


5/17/95 


10:27:41 


u 


44.00 


5. 


0 . 


879. 


0 . 


M5L17110 . IMG 


5/17/95 


10:28:01 


u 


48.60 


10. 


0 . 


0 . 


9649. 


M5L17111 . IMG 


5/17/95 


10:28:34 


h 


48.40 


2. 


0 . 


1023. 


0 . 


M5L17112 . IMG 


5/17/95 


10:28:45 


V 


48.40 


2. 


0 . 


1278. 


0 . 



103 



M5L17113. IMG 


5/17/95 


10:29:02 


u 


44.90 


5. 


M5L17114. IMG 


5/17/95 


10:29:17 


u 


48.80 


10. 


M5L17 115 . IMG 


5/17/95 


10:35:34 


h 


50.50 


10. 


M5L17 116 . IMG 


5/17/95 


10:36:05 


h 


49.60 


2. 


M5L17117 . IMG 


5/17/95 


10:36:17 


V 


49.60 


2. 


M5L17 118 . IMG 


5/17/95 


10:36:50 


u 


45.70 


5. 


M5L17119 . IMG 


5/17/95 


10:37:14 


u 


49.30 


10. 


M5L17 120 . IMG 


5/17/95 


10:38:46 


h 


49.10 


2. 


M5L17121. IMG 


5/17/95 


10:38:56 


V 


49.10 


2. 


M5L17122 . IMG 


5/17/95 


10:39:17 


u 


45.20 


5. 


M5L17123 . IMG 


5/17/95 


10:39:43 


u 


49.10 


10. 


M5L17124 . IMG 


5/17/95 


10:40:20 


h 


49.10 


2. 


M5L17125 . IMG 


5/17/95 


10:40:40 


V 


49.60 


2. 


M5L17126. IMG 


5/17/95 


10:40:59 


u 


45.70 


5. 


M5L17127 . IMG 


5/17/95 


10:41:41 


u 


45.20 


10. 


M5L17128 . IMG 


5/17/95 


10:44:21 


h 


49.30 


2. 


M5L17129 . IMG 


5/17/95 


10:44:27 


V 


49.30 


2. 


M5L17130. IMG 


5/17/95 


10:44:37 


u 


47.40 


10. 


M5L17131 . IMG 


5/17/95 


10:44:49 


u 


49.30 


10. 


M5L17132 . IMG 


5/17/95 


10:47:02 


h 


49.60 


2. 


M5L17133 . IMG 


5/17/95 


10:47:09 


V 


49.60 


2. 


M5L17134 . IMG 


5/17/95 


10:47:25 


u 


46.20 


5. 


M5L17135. IMG 


5/17/95 


10:47:39 


u 


45.20 


10. 


M5L17136 . IMG 


5/17/95 


10:47:52 


h 


49.30 


2. 


M5L17137 . IMG 


5/17/95 


10:47:57 


V 


49.30 


2. 


M5L17138 . IMG 


5/17/95 


10:48:07 


u 


45.90 


5. 


M5L17139 . IMG 


5/17/95 


10:48:16 


u 


45.20 


10. 


M5L17 140 . IMG 


5/17/95 


10:50:53 


h 


49.30 


2. 


M5L17141 . IMG 


5/17/95 


10:50:59 


V 


49.30 


2. 


M5L17142 . IMG 


5/17/95 


10:51:14 


u 


45.90 


5. 


M5L17143 . IMG 


5/17/95 


10:51:25 


u 


45.20 


10. 


M5L17144 . IMG 


5/17/95 


10:55:29 


h 


49.60 


2. 


M5L17145 . IMG 


5/17/95 


10:55:37 


V 


49.60 


2. 


M5L17146. IMG 


5/17/95 


10:55:50 


u 


45.20 


5. 


M5L17 147 . IMG 


5/17/95 


10:55:58 


u 


45.20 


10. 


M5L17148 . IMG 


5/17/95 


10:57:46 


h 


49.60 


2. 


M5L17149. IMG 


5/17/95 


10:58:06 


V 


50.10 


2. 


M5L17150. IMG 


5/17/95 


10:58:19 


u 


45.40 


5. 


M5L17151. IMG 


5/17/95 


10:58:28 


u 


45.20 


10. 


M5L17152 . IMG 


5/17/95 


11:00:33 


h 


49.80 


2. 


M5L17153. IMG 


5/17/95 


11:00:46 


V 


50.30 


2. 


M5L17154. IMG 


5/17/95 


11:01:00 


u 


45.90 


5. 


M5L17155. IMG 


5/17/95 


11:01:08 


u 


45.20 


10. 


M5L17156. IMG 


5/17/95 


11:02:55 


h 


49.80 


2. 


M5L17157 . IMG 


5/17/95 


11:03:05 


V 


49.80 


2. 


M5L17158 . IMG 


5/17/95 


11:03:20 


u 


45.90 


5. 


M5L17159 . IMG 


5/17/95 


11:03:39 


u 


44.70 


10. 


M5L17160 . IMG 


5/17/95 


11:04:22 


h 


50.30 


2. 


M5L17 161 . IMG 


5/17/95 


11:04:36 


V 


50.30 


2. 


M5L17162. IMG 


5/17/95 


11:04:58 


u 


45.40 


5. 


M5L17163. IMG 


5/17/95 


11:05:11 


u 


45.20 


10. 
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APPENDIX F. IDL IMAGE PROCESSING CODES 



The following programs are used to determine TBRD 
automatically from a file which contains all the specific 
information about a target. The desired information is 
shown in the file compilation of image names in Appendix E. 

Also the IDL code to subtract the two split-images is 
contained in this appendix. The file headers denote the use 
of the code and the code itself is commented to provide its 
own detailed explanation. 
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;fliral2.pro processes a single ship for contrast using sreadpol to get 
filenames from 
; experiment flirl2 

;Pete version that uses new contrast with no denominator 

; taken from flira8 which selects files directly 



; 2-17-95 revised xO,ix for xO always to be peak, changed shipbox to find 
peak if the 

;new peak selected with the cursor is eliminated by the threshold 

/xO,yO peak value in contiguous area 

/ix,iy, nx,ny box coordinates 

; ship temperature ship 

; shipb raw data ship 

;yc horizon line point 

;hz iy-yc location of divisor in ship box 

• ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★■a-** 

; routines readflir, shipbox, thresh, remake, fileonly, sreadpol used 
;readflir reads file data from flir format 
; horizon allows choice of horizon, formerly profiles 
; shipbox finds contiguous ship and makes ship box 

;threshld finds and returns edge threshold thresh from yy and t as well 
as peak xO,yO 

/remake uses threshold threshO to make ship and shipb from t and yy 

/fileonly separates filename from file and path for output 

/sreadpol reads polarization file and separates h or v in variable pol 



!order=0 /get rid of old program 

inverision 

iscale=2 /display 2x size 



on_ioerror, problem /if error, quit program 

.★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★•AT* 

★★★★★★★★★★★★★★★★★★★★ 



/select the files with wild cards 

pathl= ' 1 

file=pathl 

if (n_elements (path) ne 0) then print ,’ oldpath= *,path 
read, T input path for files or cr for oldpath 1 , pathl 
if(pathl ne ’’) then path=pathl+ * \ * else path=* . ’ 

if (strmid (path, strlen (path) -1, 1) ne ’ \ ' ) then path=path+ ' \ ' /fix 

oldpath error 

read, ’input filename of polarization inf ormation *, file 

sreadpol, file, files, pol /get files from flirllp 

made separately 

files=path+files /put files in standard 

findfile format 

filelist, files, path /program to list 

filenames that are different 

/for i=0,num-l do print, files (i ), i /print the filenames with 

a number 
i=* ’ 



RESTART: read, ’ Input start and stop number for processing CR to stop*,i 

if strlen (i) eq 0 then goto, finish 
istart=fix (i) 
istop=istart 

read, 'Stop number 1 , istop 
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C5N=DBLARR ( (ISTOP-ISTART+1 ) , 3 ) ; The first z is full box, second, above 

horiz, third below 

read, 1 input first guess difference t level, eg 1 to include all of 
ship 1 , tdif 

r 

★ ★★★★★★★★★★★★★★★★-A-** 

close, 2 

fileout=* * 

read, ' input output filename for contrast results ', fileout 
openw, 2, fileout 

printf , 2 , systime ( ) , * date of analysis* 
printf , 2, istop-istart+1, ' number of files* 

r 



auto=* * 

read, * input CR for semi-automatic proccessing * , auto 

tdif save=tdif ; save tdif so that changes will 

not affect batch 

printf ,2, ‘Temperature difference = *,tdif ; save in a file 

; tdif, nx , ny , yc , iy 

printf , 2, * filel, tdif, nx, ny, yc (horizon) , iy (ship bottom) , Polarization* 
for i=istart, istop do begin ;ii = first of pair 

tdif=tdifsave /restore the temperature to the 

original batch value 
f=files (i) 



file=f ; save the name 

readf lir, f, yy, t, bar, tt, head /read flir data, convert to 

temperature, add bar 

• *************************** * * clisplay the raw 
da t a ************************** ** ***** ★ *★★★*★ * 



loadct,0 /set color table 0 

iscale=2 /make pictures twice their normal 

size 



; remove the path length from the file string and label the tv 
image with the filename 

pathlen=strlen (path) +1 
filelen=strlen (file) 
labelen=filelen-pathlen 
label=strmid ( file, pathlen, labelen) 

window, 0, xsize=140*i scale, ysize=140*iscale+50, xpos=0, ypos=0, title= * 0 
' +label 

tvscl, rebin (tt, 140*iscale, 165*iscale) /display the image 
• ★ End of temperature 

sp lay * 

window, 2, xpos=0, ypos=165*iscale+35, xsize=140*iscale, ysize=140*iscale 



print, * choose horizon with left button or reject picture with right* 
tvscl, rebin (yy, 140*iscale, 140*iscale) 
binary image raw 

HORIZON, smooth (rebin (yy, 140*iscale, 140*iscale) , 5 ) , xc, yc 
position 

if(!err eq 4) then goto, skiploop 
processing of file 
xc=xc/iscale 
yc=yc/ iscale 
size not picture size 
print, *xc,yc= *,xc,yc 

; ****************************thresh 0 ld with temperature 
f f g^encs*************************** 



/display the 
; read cursor 
/skip 

/scale to 140 
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threshld, t, yy, thresO , xO, yO 
thres and find xQ,yO peak 



; define 



; RESHIP RESTORES THE VARIABLES MODIFED BY SHIPBOX, NOT THRESH 
errorct=0 ; don’t allow infinite 

loopes of choices 

• ★ RE SHI p**************** - *''*''*''*''*''*'*''*''*''*' 



reship: 

remake, ship, shipb, t, yy, thresO 
thresholded ship and shipb using thresO 
tvscl, rebin (shipb, 140*iscale, 140*iscale) 
image raw 

tvcrs, xO*iscale, yO*iscale 

shipbox, shipb, ship, xO, yO, ix, iy,nx,ny, i scale 
ship box contiguous 



; remake the 

; display the binary 

; cursor on hot spot 
;ix,iy,nx,ny define 



; WARNING ALL VARIABLES EXCEPT I SCALE ARE CHANGED IN THIS 



CALL 

boxdraw, ix*iscale,iy*iscale,nx*iscale,ny*iscale ; show where the box 
is in *iscale 

tvcrs, xO*iscale, yO*iscale, /device ;put the cursor on 

the peak location 

• ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★TIT 
r 

'k'k-k'k’k'k'k'k'k'k'k'k'k'k'k'k'k-k'k'k'k 



; *** a k ove makes box with minimum 

thrsshold*************************************************** 

! ERR=4 

; if (auto ne ’ ’ or errorct gt 0) then begin ; suspend auto only 

for errors 

print, ’right button, use auto position, left, new cursor position* 
cursor, ixnew, iynew, /device, /down ;when pushed, 

choose 

print, * !ERR= ’ , ! ERR 

if ( ! ERR eq 4) then print, ’auto selected* 
if(!ERR eq 1) then begin 
xO=ixnew/iscale 
y0=iynew/i scale 

print, ’new peak at ’,x0,y0 ;NEW peak 

at x0,y0 

remake, ship, shipb, t, yy, thresO ;if not 

remade, the shipbox will have destroyed shipb 

/WARNING SHIPBOX DESTROYS SHIPB AND SHIP 

shipbox, shipb, ship, xO, yO, ix, iy, nx, ny, i scale ; ix, iy,nx, ny 

define ship box contiguous 

boxdraw, ix*iscale, iy*iscale, nx*iscale, ny*iscale /show where 

the box is in *iscale 
endif 

/endif / semi auto is 

two button pushes now 



/ ***************** IF THE hot SPOT IS FOUND DO LOCAL THRESHOLD IN THE 

remake, ship, shipb, t, yy, thresO+tdif / remake with full 

threshold 

print, ’ redisplay with full threshold* 

shipbox, shipb, ship, xO, yO, ix, iy, nx, ny, i scale / xO, yO, nx, ny 

define ship box contiguous 

tvscl, rebin (shipb, 1 40 *i scale, 1 40 *i scale) 

boxdraw, ix*iscale, iy*iscale, nx*iscale, ny*iscale 
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;place the cursor on 



tvcrs, ix, iy 
the max for a check 
print, ’cursor on maximum* 

. ************************* fi n d a new shipbox based on threshold and 



if (total (shipb) eq 0) then begin ;if data is all zero, 

threshold did not work 

errorct = errorct+1 /avoid an infinite 

loop 

print, 'no ship was found. Data is not usable with this tdif',tdif 
if (errorct gt 2) then begin 
print,' skipping ' , file 
goto, skiploop 
endif 

read, ' input new tdif',tdif 

; auto= ' ' /suspend auto 

permanently on first error 
goto, reship 



endif 



tdif f= 1 * 

print, ' tdif= ' , tdif 
tdif f= ' ' 

if (auto ne '' or errorct gt 0) then $ 

read, 'input new tdif or CR for old or negative to skip', tdif f 
if (strlen ( tdif f ) eq 0) then tdiff=tdif else begin 
tdif=float (tdif f ) 

if (tdif It 0) then goto, skiploop 
goto, reship 
endelse 



• ★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★★ 
r 

★★★★★★★★★★★★★★★★★★★★ 



/choices section checks to see which picture is larger and uses the 
appropriate tdif 

/to remake and replot the temperatures reference windows 

/ ship temperature ship in box & sv extension 
/box temperature box not ship 
/ shipb binary ship in box 
/boxb binary box 



ship=(t gt thresO + tdif ) *t 
the OLD threshold 

shipb=(t gt thresO + tdif ) *yy 



/be sure NOR to use 



ship=ship (ix: ix+nx, iy: iy+ny) 
shipb=shipb (ix : ix+nx, iy : iy+ny) 
box=t (ix: ix+nx, iy: iy+ny) 
box from original data 

boxb=yy (ix : ix+nx, iy: iy+ny) 
box=box* (ship eq 0) 
(temperature) in box 



/ subset 

/everything in the 



/not ship 
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boxb=boxb* ( shipb eq 0) 
zscale=140/nx 

window, 5, xpos=140*iscale, ypos=0, xsize=nx*z scale, ysize^ny* z scale 
tvscl , congrid (ship, nx*zscale, ny* z scale) 

******************** 

; make contrasts ell = first file, shipfull, cl2=above horizon, cl3=up 
to horizon 

; C3N full box, above horizon, below horizon 
;yc horizon line point 

; iy ship bottom 

; hz yc-iy location of divisor in ship box 

★★★★★★★★★★★★★★★★★★★ 

hz=yc-iy 

;if(hz It 0) then horizon is below ship and ship is flying 
;if(hz ge iy+ny) then horizon is completely above ship 

c5n (i-istart, 0) =pcontra (shipb, boxb) 

; f rom yc (horizon) -iy to ma 

if (hz le 0 or hz ge ny) then begin ;hz 0 means c5n(i,2) 

is ill defined as ship is on horizon 

print, 'hz,ny f ,hz,ny, 1 error' ;hz gt ny means 

c5n(i,l) is below horizon 

c5n (i-istart, 1) =c5n (i-istart , 0) ;hz It 0 means ship 

is flying above horizon 

c5n (i-istart , 2) =c5n (i-istart, 0) 
endif else begin 

c5n (i-istart, 1 ) =pcontra (shipb (*, hz : *) , boxb ( * , hz : * ) ) 
c5n (i-istart, 2) =pcontra (shipb (*, 0:hz) ,boxb (*, 0:hz) ) 

; duplicate the first if the value hz is zero 



endelse 

print, fileonly ( files ( I) , path) , * ' , tdif , nx, ny, yc, iy, ' ',pol(i) 

print, C5N( (I-ISTART) ,0) , f ormat=" ( (F7.3,2x) , 'Full ship' ) " 
print, C5N( (I-ISTART) , 1) , format=" ( (F7. 3,2x) , 'Above Horizon' ) 11 
print, C5N( (I-ISTART) ,2) , format=" ( (F7. 3, 2x) , 'Below Horizon' ) " 

print f , 2, fileonly ( files (I ) , path) , ' ' , tdif , nx, ny, yc, iy, ’ * ,pol (i ) 

printf , 2, C5N ( (I-ISTART) ,0) , format^" ( (F7. 3,2x) , 'Full ship' ) " 
print f, 2, C5N ( (I-ISTART) ,1) , format='' ( (F7.3,2x) , 'Above Horizon' ) " 
printf, 2, C5N ( (I-ISTART) ,2) , format=" ( (F7. 3,2x) , 'Below Horizon' ) " 



skiploop : 
endfor 
/ 

★*★*★*★***★★*★*★**★* 
goto, finish 

problem: print, 'I had a problem and quit' 
finish: print, 'program finished' 
close, 2 
end 

/program shipszl2 reads output of fliral2 (Pete Heisey) which is: 
/filename, tdif ,ix, iy,nx,ny,H or V 



130 



/and calculates TBRD with ship size and plots it 

close, 1 ; shut up in case of an error 

tempp= 1 



file=' 1 

read, 'input filename from analysis output', file 

on___ioerror, stoperr 

openr, 1, file 

nfiles=0 

iy=0 

yc=iy 

nx=iy 

ny=iy 

i=0 

t=l . 0 



HV= ' 



;file count 
; read to end of 



file 



•read string and print 



;if number of files not 



while not eof(l) do begin 
which is not more than nfiles 
; loop : readf , 1 , tempp & print, tempp 
it 

loop : readf, 1 , tempp 
if (nfiles eq 0) then begin 
fount yet, look for it 

if ( st rpos (tempp, ' number of files') gt 0) then begin 

reads, tempp, nfiles & print , tempp, nfiles ;get the number of files 

rdata=fltarr (5, nfiles ) ;make array of number of 

files, may not use it all 
cdata=fltarr (nfiles, 3) 

; namearr=strarr ( 12, nfiles) 
namearr=strarr (nfiles ) 
pol=strarr (nfiles) 

print, ' finished making arrays, Nfiles= 
know arrays were made 
endif 
endif 

if (strpos (tempp, IMG' ) It 0) then goto, loop else $ 

namearr (i ) = (strmid (tempp, 0, strpos (tempp, ' G* ) +1 ) ) /search for file end 
IMG 



/array of contrast data 
/array of filenames 



, nfiles 



/let me 



/print, (strmid (tempp, strpos (tempp, ' G' ) +1, strlen (tempp) ) ) 

y=strcompress (strmid (tempp, strpos (tempp, ' G' ) +1, strlen (tempp) ) ) / get 

rid of white spaces 

/ reads, y, t, ix, iy, nx, ny, HV & print , t, ix, iy, nx, ny, HV /convert data and 

print it 

/printf , 2, tdif , nx, ny, yc, iy maee this 

reads, y, t , nx, ny, yc, iy, HV 

rdata (0, i) =t 

rdata ( 1, i ) =nx 

rdata (2 , i ) =ny 

rdata ( 3, i ) =yc 

rdata (4, i ) =iy 

pol (i)=HV 

y=l . 0 

readf, l,y & cdata(i,0) = y 

readf, l,y & cdata(i,l) = y 

readf, l,y & cdata(i,2) = y 

print ,' incrementing loop counter' 
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; increment the stored 



i=i+l 
variable 

if ( rdata ( 4, i-1 ) le 1) then i=i-l ;if height is 0, skip 

the data 
endwhile 
goto, stopp 

stoperr: print, 1 io error 1 

goto, loop 

stopp: close, 1 

if(i It 0) then i=0 

print, 1 finished 1 , i 

rdata=rdata (*, 0 : i-1) 

cdata=cdata ( 0 : i-1, * ) 

namearr=namearr (0 : i-1) 

print, 1 file= 1 , file 

S PLOT 9e , RDATA, CDATA, FILE , pol 

end 



/program contrast. pro which calculates TBRD for each ship image 
/anayzed. it is a subroutine of shipszl2 

function contrast, ship, box 
contrast= ( shipav ( ship) -ship (box) ) 
return, contrast 
end 



/program splot 9 is subroutine of shipszl2 which takes the ship 
conatrast data and plots it vs. ship range determined by size 

PRO SPL0T9, RDATA, CDATA, file, pol 
ptitle=file 

;cdata(*,0) contrast for whole ship 

;cdata(*,l) contrast for above horizon 

/cdata(*,2) contrast for below horizon 

/rdata (2,*) ship size from height 

hpol=where ( strtrim (pol, 2) eq ? h r ) /array of number for h 

polarization 

vpol=where ( strtrim (pol, 2 ) 
help, hpol , vpol 
if(n elements (hpol) le 1) 



eq ? v ! ) /array for v polarization 

then hpol=where (strtrim (pol, 2) ne 'h 1 ) 



132 



if (n_elements (vpol) le 1) then vpol=where (strtrim (pol, 2 ) ne 'v') 

LOADCT, 39 

;read, 'Input title ', ptitle 
; WINDOW, 0,XSIZE=1024, YSIZE=768 

; plot, 1. /rdata (2, *) , cdata (*, 2 ) , CHARSIZE=2 , psym=2 , color=0, background=255, 
charthick=2, title=* ! 17 * +ptitle $ 

; , xtitle= ' Ship Range ', ytitle= 1 Target-Background Radiance 
Difference ' , /nodata 

; oplot, 1 . /rdata (2, *) , cdata (*, 2 ) , psym=2, color=0 
; oplot, 1 . /rdata (2 , * ) , cdata ( * , 0 ) , psym=l, color=0 
; oplot, 1 . /rdata (2 , * ) , cdata (*, 1) ,psym=4, color=60 
;xyouts, .65, .87, charsize=3, charthick=2, ' ! 17Above 
Horizon 1 , color=60, /normal 

;xyouts, .65, .82, charsize=3, charthick=2, * ! 17 Full Ship* , color=0, /normal 
;xyouts, .65, .77, charsize=3, charthick=2, ' ! 17Below 
Horizon ' , color=60, /normal 
; WINDOW, 1 , XSI ZE=1024, YSIZE=7 68 

;plot, 1. /rdata (2, *) , cdata (*, 0) , CHARSIZE=2 , psym=2 , color=0,background=255, 
charthick=2, title= f ! 17 ' +ptitle $ 

; , xtitle=' Ship Range 1 , ytitle= 1 Target-Background Radiance Difference' 

; oplot, 1. /rdata (2, *) , cdata (*,1) , psym=4, color=0 
; ; oplot, 1 . /rdata (2 , * ) , cdata ( *, 2 ) , psym=6, color=0 
WINDOW, 2 , XSI ZE=102 4 , YSI ZE=7 68 

plot, 1 . /rdata (2, * ) , cdata (*, 2 ) , CHARSIZE=2, psym=2 , color=0, background=255, c 
harthick=2, title= '117' tptitle, /nodata $ 

, xtitle= 1 Ship Range from Ship Size ', ytitle =t Target-Background Radiance 
Difference 1 

oplot, 1 . /rdata (2 , hpol) , cdata (hpol, 2 ) , psym=2 , color=0 
oplot, 1 . /rdata (2 , vpol ) , cdata (vpol, 2 ) , psym=6, color=0 
oplot, 1 . /rdata (2 , hpol) , cdata (hpol, 1 ) , psym=l, color=0 
oplot, 1 . /rdata (2, vpol ) , cdata (vpol, 1 ) ,psym=4, color=0 
xyouts, . 65, . 77, charsize=l, charthick=2, 1 ! 17 *+ = H ! 9B ! 9V ! 17= 

V* , color=0, /normal 

xyouts, . 65, . 87 , charsize=l, charthick=2, 1 ! 17Above Horizon * 

! 9B* ,color=0, /normal 

xyouts, . 65, . 82, charsize=l, charthick=2, '! 17Below Horizon + 

! 9V 1 , color=0 , /normal 
end 



;this program subtracts the arrays of the two split images from the 
splitimage polar ime ter. 

/program splitcom compares split image in yy array 

getfiles, files 

/get array of filenames in files 

loop: 

/do it again 

read, 'input number of file to get f ,ifile 
if(ifile It 0) then goto, problem 
f=files (ifile) 

readflir, f , yy, t, bar, tt, head ; get 

file and make temperature array 

! order=Q 

on_error, 1 

ybot=yy ( 0 : * , 7 0 : * ) 

ytop=yy ( 0 : * , 0 : 69 ) 
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; Invert the image 



yes= 1 1 

read, ’ input y to invert’, yes 
if (yes eq ’y') then begin 
ybot=ybot (*, 69-indgen (70) ) 
ytop=ytop (*, 69-indgen (70) ) 
endif 

hotspot, ytop, xhott, yhott ; FIND THE 

HOTSPOT IN THE IMAGE 



hotspot , ybot, xhotb, yhotb ; subroutine 

window, 1, xsize=140*3, ysize=140*3 
tvscl, rebin (ybot, 140*3,70*3) 

tvscl, rebin (ytop, 140*3, 70*3) , 0, 70*3 ;display both images 

xdif=xhott-xhotb 

; calculate displacement to overlay hotspots 

ydi f =yhott-yhotb 

print, xdif, ydif, * xdif ydif’ 

ytop=float (ytop ) ;make 

floating so subtraction will work 
ybot=float (ybot) 

displace, ybot, fix (xdif ), fix (ydif ) ; subroutine which 

; overlaps images 

window, 3, xsize=140*4, ysize=70*4 

absoldif=abs (ytop-ybot) ;make the 

values of the difference positive 

tvscl, rebin (absoldif , 140*4,70*4) ; display difference 

goto, loop 

problem: print, ’exiting 1 

; error and normal exit for -1 filename 



end 



;pro displace, yval, xdif , ydif 

ysize=size (yval) 

print , ysize 
if (ydif ge 0) then 

yval ( 0 : * , abs (ydif) : ysize (2 ) -1 ) =yval (0 : * , 0 : ysize (2) -1-abs ( ydif) ) $ 
else yval (0: *, 0 : ysize (2 ) -1-abs (ydif) ) =yval (0 : *, abs (ydif) : ysize (2 ) -1 ) 
end 

pro hotspot , yval ,xhot,yhot 

ysize=intarr (5) 
ysize=size (yval) 

print , max (yval, ymax) & print, ymax 

print, ymax- (ymax/ysize (1) ) *ysize (1) , ymax/ysize (1) 

xhot= ymax- (ymax/ysize (1 ) ) *ysize ( 1 ) 
yhot=ymax/ysize ( 1 ) 
end 
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APPENDIX G. 



SPLIT IMAGE POLARIMETER IMAGES, 3 AUGUST 1995 



The following appendix contains the images recorded of 
the USS John McCain (Arleigh Burke class DDG) on 3 August, 
1995 with the split-field polarimeter. This was a 
preliminary experiment conducted to test the performance of 
the optical system in the far field and to validate field 
performance of all system components. Any Data on the image 
file not contained in the Experimental notes is stored in 
the image file header. 
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Experimental Notes 



Image File 
Name 


Targe 

t 


Thermal 

Range 


Remarks 


m0803001 


BB 




setup test-BB=black body source 


m.0 8 030 02 


BB 




infinite focus, BB @20' 


m0803003 


BB 




ff 


m0803004 


BB 




Black body temperature 36degrees 


m0803005 


McCain 


2 


cal table 4003 used 


m080300 6 


McCain 


5 




m0803007 


McCain 


2 




m0803008 


McCain 


2 




m0803009 


McCain 


2 




m0803010 


McCain 


5 




m0803011 


McCain 


2 




mO 8 03 01 2 


McCain 


5 




m0803013 


BB 


2 


cal tab. 4004, BB 36 degrees 
OC=-44 


m0803014 


BB 


5 




m0803015 


McCain 


5 




m080301 6 


McCain 


2 




m0803017 


McCain 


2 




m0803018 


McCain 


2 




m0803019 


FOV 


2 


empty FOV with aft of ship on 
far right side 


m0803020 


McCain 


5 
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m0803021 


McCain 


10 




m0803022 


McCain 


5 




m0803023 


BB 


10 


focus set on inf. for BB 
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I 
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